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ABSTRACT 


Polycystic ovarian syndrome (PCOS) is the most prevailing endocrine 
disorder among women of reproductive age, worldwide. Its prevalence is 
52% in Pakistani women and the rate of infertility is gradually increasing. 
It can be treated through early diagnosis, management of symptoms, and a 
healthy lifestyle. This cross-sectional study aimed to determine the level of 
awareness and knowledge regarding PCOS in the women of Punjab, 
Pakistan. For this purpose, a total of 350 women between the ages of 16 and 
50 years were included in the current study. The data was obtained through 
a questionnaire and analyzed through SPSS (version 25). The results 
determined that 5696 of women were found to have the knowledge of PCOS, 
while 44% had never heard about it. Furthermore, about 67% of 196 PCOS- 
aware women were unaware of its signs and symptoms. Moreover, in 1496 
of the cases, a friend was the source of knowledge regarding PCOS. 
Whereas, approximately 10% of these women had heard about the 
syndrome through education, 9% through their family, 7% through various 
media such as newspapers, television, and social media, and 17% through 
other sources including general physicians and gynecologists. It was 
concluded that only 33% of the participants who claimed to have heard of 
PCOS were aware of its symptoms, whereas 67% were unaware of them. A 
total of 152 women, out of 350, had never heard about PCOS. 


Keywords: endocrine disorder, infertility, healthy lifestyle, polycystic 
ovarian syndrome (PCOS), women of reproductive age 


1. INTRODUCTION 


Polycystic ovary syndrome (PCOS) is known worldwide as a common 
syndrome affecting the women of reproductive age and causing 
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reproductive, endocrine, and metabolic abnormalities. PCOS is delineated 
as a heterogeneous disorder associated with hormonal changes, a deranged 
reproductive system, and complications related to the metabolic system [1]. 
Itis a heritable disorder which affects women throughout their lifetime and 
can be transferred to the next generation, although the underlying genetic 
factors remain unknown [2]. Hyperandrogenism, ovaries with multiple 
cysts, and impaired ovaries are the key features of PCOS. This disorder 
affects both developed as well as developing countries. Its prevalence 
ranges from 4-12% in the general population of women of reproductive age 
across the world [3]. 


Throughout puberty, the prevalence of PCOS continuously increases 
and reaches 26% by the age of 15. Insulin resistance or pre-diabetes and 
Type 2 diabetes are prevalent in women with PCOS living in Asia. Insulin 
resistance and high levels of insulin are two common precursors of Type 2 
diabetes and have a high prevalence in South Asian women. Type 2 diabetes 
is inherited in the form of a complex genetic trait that interacts with other 
factors and contributes to developing PCOS [4]. Symptoms commonly 
associated with PCOS are irregular or missed periods due to anovulation, 
ovaries consisting of multiple cysts, non-appearance of ovulation, high level 
of androgens [5], oligomenorrhea, amenorrhea, prolonged erratic menstrual 
bleeding [6], infertility [7], hirsutism [8], acne [9], mood swings including 
depression and anxiety [10], and family history with Type 2 diabetes [11]. 


Three distinct clinical phenotypes were suggested by the Rotterdam and 
AE-PCOS Society namely Frank PCOS (oligomenorrhea, higher level of 
androgens, and PCO), ovulatory PCOS (PCO, higher level of 
hyperandrogenism and normal menstruation), and non-PCO PCOS 
(oligomenorrhea, hyperandrogenism, and normal ovary) A fourth 
phenotype, that is, mild or normoandrogenic PCOS is also recognized by 
the Rotterdam criteria, which is characterized by oligomenorrhea, PCO, and 
typical androgen concentrations [12]. 


In addition to other classical symptoms, stress is highly coexistent with 
PCOS. Stress in PCOS can be of various types namely psychological 
distress, mental disorders, anxiety, emotional disarray, and depression [10]. 
Along with definitive mental stress, ‘metabolic stress’ is also quite common 
in women with PCOS. Indeed, metabolic stress is the most influencing 
stress in PCOS. It is often referred to as the stimulant behind all the 
emotional and anxiety-related consequences in women with PCOS [13]. 
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A characteristic ofthe ovaries in PCOS is rapid transition from primitive 
follicles to growing follicles with rising numbers of follicles from 2 mm to 
3 mm and then from 3 mm to 4mm [14-16]. Adolescent girls with PCOS 
features often require the management of their symptoms before the 
diagnosis of PCOS. The management of adolescents with a confirmed 
diagnosis of PCOS should include its awareness and lifestyle interventions 
[17]. The diagnosis of girls with signs and symptoms indicating PCOS 
includes family history and a complete physical examination. Typically, 
laboratory examination includes thyroid profile as well as the determination 
of prolactin, total testosterone, androstenedione, sex hormone-binding 
globulin (SHBG), and dehydroepiandrosterone (DHEA). Furthermore, 
adrenal imaging tests can be used for diagnosis depending on clinical 
examination, physical examination, and initial laboratory results [18]. 


Mohamed conducted a quasi-experimental study on 300 women that 
included implementing and designing a helpful instructional guideline 
about endometriosis. The study shows an important result, indicating an 
improvement in females' awareness about it. [19]. Belenkaia (2019) showed 
that polycystic ovarian syndrome (PCOS) is recognized as one of the most 
common endocrine diseases in women, universally. On the other hand, this 
term does not capture the variety of medical symbols linked with this 
syndrome. Menstrual abnormality and medical features of androgen excess, 
which are common place in women with PCOS, are not included under the 
definition of PCOS, which is limited to polycystic ovarian morphology 
(PCOM) [20]. 


2. METHODOLOGY 
2.1. Study Design 


A cross-sectional study was conducted focusing on assessment-based 
questions and the provision of education about PCOS. 


2.2. Locality Setting 
Data was collected from the various cities of Punjab, Pakistan. 
2.3. Duration 
This study was completed in October 2021. 
2.4. Sample Size 
A total of 350 questionnaires were distributed among women. 
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2.5. Sample Technique 
It is a questionnaire-based study. 
2.6. Inclusion Criteria 


Women between the ages of 16 and 50 years were included in the 
current study. 


2.7. Exclusion Criteria 
Women aged less than 16 years or more than 50 years were excluded. 
2.8. Data Collection Procedure 


Data was collected through a questionnaire developed beforehand. The 
questionnaire was constructed in parts related to demographics, information 
sources, understanding of signs and symptoms, risk factors, prevalence, and 
etiology. 


2.9. Data Analysis 


Data was tabulated and analyzed through Microsoft Excel and IBM 
SPSS Statistics (version 25). Chi-square test was used to statistically 
analyze the comparison between categorical variables, that 1s, variables that 
only take a finite number of distinct values. The results were reported in 
terms of frequencies and percentages calculated for each question in the 
questionnaire. 


3. RESULTS 
3.1. Demographic Characteristics 


Table 1. Demographic Characteristics 


Description Frequency Percentage 
Age 
15-25 years 196 56 
26-35 years 78 22 
36-45 years 41 11.7 
46-55 years 35 10 
Locality 
Urban 250 28.6 
Rural 100 71.4 
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Description Frequency Percentage 
Education Level 
School 45 12.9 
College 162 46.3 
University 141 40.3 
Other 2 6 
Marital Status 
Married 152 43.5 
Single 198 56.6 
Occupation 
Employee T] 22,5 
Household 90 25.5 
Student 183 52.3 


Demographic characteristics are tabulated in Table 1. The table shows 
that 196 (56%) women were from the age group of 15-25 years, 78 (22%) 
were from the age group of 26-35 years, and 41 (11.7%) were from the age 
group of 36-45 years. It also shows that 100 (28.6%) women belonged to 
the rural areas of Punjab, while the rest of the 250 (71.4%) women belonged 
to urban areas. Furthermore, 162 (46.3%) women were enrolled in college 
and 141 (40.3%) were university going. Whereas, 198 (56.6%) were single. 


3.2. Knowledge Regarding PCOS 


Knowledge regarding PCOS is tabulated in Table 2. The majority of the 
women, that is, 56% responded positively to questions regarding PCOS 
awareness. However, 44% had never heard of PCOS and did not know 
anything about it. 


3.3. Diagnosis of PCOS 


The number of women diagnosed with PCOS is tabulated in Table 3. It 
shows that 82.9% of women had not been diagnosed with PCOS, while 
17.1% of them had been diagnosed with the disease. 


The results in Table 4 show that there is a significant association 
between PCOS diagnosis and menstrual problems in women (p-value is less 
than 0.05 which is significant). According to the above table, 32 women 
faced skipped periods problem and 16 faced irregular menstrual cycles. 
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Table 2. Knowledge and Awareness regarding PCOS 


: Do not 

Questions Yes No Enn 
Have you ever heard of PCOS? ( f ) ( ns ) N/A 
Do you know the symptoms of 116 234 N/A 
PCOS? (33.1%) (66.994) 
Do you think PCOS is linked with 70 81 199 
Diabetes Mellitus Type 2? (20%) | (2394) (56.9%) 
Do you think PCOS has a negative 
impact on the quality of life in the on os 
form of stress, anxiety etc? Vegan) (uu) 027) 
Do you think untreated PCOS can 94 59 197 
cause infertility in women? (26.9%) (16.9%) (56.3) 
Do you think obesity increases the 119 54 177 
risk of developing PCOS? (34%) (15.4%) (50.6%) 
Do you think unpleasant mood is 105 62 183 
also a risk factor for PCOS? (30%) (17.7%) (52.3%) 
Is ultrasound a best way for the 97 59 194 
detection of PCOS? (27.7%) (16.9%) (55.4%) 

Table 3. Diagnosis of PCOS 
Description Yes No 

Have you ever been diagnosed with PCOS? 82.994 17.1% 


Table 4. Association between PCOS Diagnosis and Menstrual Problems in 


Women 
History of PCOS Diagnosed Patients with Menstrual Problems 
Irregular Skipped . Clots of 
oneor Excessive : 
Yes/No menstrual i blood in p 
p more bleeding 3 
periods : periods 
periods 
Have you 
ever been No 290 0 0 0 0 
diagnosed 0.00 
with Yes 0 16 32 5 7 
PCOS? 
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Table 5. Association between WHO Obesity Criteria and Number of 
Diagnosed Patients 


Have you ever been diagnosed 


with PCOS? p-value 
Yes No 
Underweight 36 2 
WHO Obesity Normal 199 11 0.00 
Criteria Overweight 44 39 i 
Obese 11 8 


The results show that there is an association between WHO obesity 
criteria and the number of diagnosed patients (p-value is less than 0.05 
which is significant). Following the WHO obesity criteria, 36 underweight, 
199 normal, 44 overweight, and 11 obese women were diagnosed with 
PCOS. 
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18-25 26-35 36-45 46-55 
Age Groups of Women 


Figure 1. Age Groups of Women and their Education Level 


Figure 1 depicts that most women in the age group 18-25 were college 
educated. While, most women had university education in the age group 26- 
35. Whereas, most women were college educated in the age group 36-45. 
Finally, most women in the age group 45-55 were also college educated. 
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Through Through Through Through Other None of 
friends family education media sources these 


Source of Information regarding PCOS 


Figure 2. PCOS Awareness in Women and their Source of Information 


Figure 2 depicts that out of the 196 women who had heard of PCOS, 48 
heard about it through friends, 30 through family, 32 through education, 23 
through media, and 59 through other sources. 
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Have you Doyouhave Have you Do you Do you get 
been irregular skipped one or experience clots of blood 
diagnosed menstrual more periods? excessive in your 
with PCOS periods? bleeding? periods 


History of PCOS Diagnosed Patients with Menstrual Problems 


Figure 3. PCOS Diagnosis and Menstrual Problems in Women 
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Figure 3 depicts that 32 women had faced skipped periods problems and 
16 had faced irregular menstrual cycles. 


4. DISCUSSION 


Sunanda recorded that 8596 of women (in their study?) were in the age 
group of 21-25 years, while 82% of them enjoyed a diet incorporating 
diverse food groups, and 92% had an orderly menstrual pattern [21]. Sills 
noticed that more than half (6396) of 657 women involved were in the age 
range of 26-34 years at the time. Approximately 60% of them were normal 
according to the WHO criteria of obesity. Nearly 2496 were overweight and 
6% fell in the obese category. Conversely, 11% of women were 
underweight [22]. Sanchez noticed that 32% of women were obese in their 
study [23]. 


Moreover, about 5696 of women under observation had heard of PCOS. 
On the contrary, 44% of them had never heard of the disease and were 
unaware of its symptoms, risk factors, treatment options, and severity of 
consequences. Furthermore, 67% of women, who had heard of PCOS, did 
not know about its symptoms. Approximately 33% of women were aware 
of the symptoms of the disease. A survey by Scott estimated that more than 
50% of the 657 women involved were fully aware of the syndrome. In this 
regard, as shown by the analysis, doctors were the primary source of 
knowledge for 63.5% of women [23]. 


Out of the 196 women who had heard of PCOS previously, the source 
was a friend in 14% of cases. Whereas, 10% of them had heard of the 
syndrome through education, 9% through their family, and 7% through 
different media including social media, newspaper, and television. 
Furthermore, 1796 of women said that their source of information was their 
general physician or gynecologist. Sunanda observed poor awareness of 
PCOS in 76% of women, whereas 10.7% had thorough information of the 
disease [21]. Sills concluded that over 97% of respondents were acquainted 
with POCS, while 1.9% had not been briefed about the disease [23]. Sills 
also noticed that all respondents who were 26 to 34 years of age were 
comparatively more informed of the condition than those belonging to other 
age groups [23]. Gul reported that 20 out of 177 participants were conscious 
of this condition. Of these 20 participants, 11 had qualifications in 
biological sciences [24]. 


In the current study, only 14% of women considered irregular menstrual 
y, only g 


10 SF R Scientific Inquiry and Review 
dag Volume 7 Issue 1, 2023 


Rafique et al. 


cycles as something which should not be neglected. Almost 51% didn't 
know that obese women are at a higher risk of developing PCOS than 
others, whereas 15% said that obesity has no role in causing the disease. 
Almost 53% of women didn't know that unpleasant mood is a risk factor 
for PCOS, while 18% said that it is not related to the disease. Approximately 
17% of women didn't know that infertility can be induced because of an 
untreated PCOS. Only 15% identified hirsutism to be a condition that 
cannot be ignored and that could lead to PCOS. These problems, that is, 
menstrual irregularity, hirsutism, and obesity, which are often neglected by 
women, are associated with elevated testosterone levels and PCO 
appearance on ultrasound [22]. Joshi concluded that a record of infrequent 
menstrual periods showed a 93.3% positive probability of the disease and 
an 86.7% negative probability to identify its potential occurrence [24]. 


Among 60 women diagnosed with PCOS, 41.7% believed in physical 
activity. Approximately 31% were aware that physical exercise can 
improve PCOS related condition. Many women with the disease experience 
the feelings of hopelessness and insecurity because of their inability to 
conceive, the size of their bodies, and excessive facial hair. This has a huge 
effect on the psychological health of women with PCOS. Scott’s research 
suggested distress in 66.5% of women due to PCOS, while this analysis 
found 32% of women as apprehensive because they were not fully informed 
about the nature of the disease and if they would be able to control it. About 
19% of women feel distressed as the out of form body and the decreased 
quality of life obstruct their emotional stability, which is also affected by 
chemical shifts in PCOS [23]. Maiya stated that ranked resistance training 
helped to reduce the cyst’s thickness, improved fertilization, as well as the 
conceiving rate, as compared to an experimental group [25]. A similar 
survey by Nidhi found that therapeutic meditation system in PCOS patients 
minimizes the feelings of depression [26]. 


In the current study, about 64% of women didn’t know that androgen 
levels are significantly raised during PCOS. Furthermore, approximately 
12% chose unaltered androgen, whereas 9% chose low androgen levels as 
blood test confirmation for the disease. The androgen abundance and PCOS 
community developed the newest testing manual in 2006. According to this 
manual, hyper-androgenism and dysfunctional ovaries (oligoavulation or 
anovulation and/or cysts in ovaries) are the salient features of the disease 
[27]. Approximately 40% of women were positive that the treatment of 
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PCOS is possible. Nearly 12.596 declared it as an untreatable disease. About 
47.5% didn't know if any treatment for PCOS exists or not. 


The treatment of PCOS can be initiated through increased physical 
activity and regular exercise. Changes in routine are regarded as the first 
line therapy for women diagnosed with the disease [28]. In a relevant study, 
Sills found that 99% of women showed interest in other possible treatments 
of their condition, rather than using fertility medication or contraceptives 
[23]. Unawareness in women regarding PCOS is a result of the lack of 
knowledge and poor education about this matter. Illiteracy is also one of the 
leading causes of the lack of awareness about such a sensitive issue. 


4.1. Conclusion 


Based on the survey findings, it was concluded that a significant portion 
of women, approximately 3396 of those who claimed to have heard of 
PCOS, demonstrated an understanding of the syndrome and its associated 
symptoms. However, it was concerning to note that the majority, 
constituting 67% of these women, lacked awareness regarding the signs and 
symptoms of PCOS. Moreover, a substantial number of participants, 
specifically 152 out of the total 350 women surveyed, acknowledged never 
having heard of PCOS before. 


Healthcare professionals should play a crucial role as primary sources 
of information about PCOS for women. Targeted educational campaigns 
can bridge the knowledge gap and enhance awareness, leading to improved 
understanding, early detection, and management of PCOS. 


4.2. Recommendations 


Exercise plans that are well-designed for PCOS patients should be 
introduced. Physicians should promote regular physical activity apart from 
providing need-based exercise plans. Gynecologists should coordinate with 
physiotherapists to ensure that PCOS patients are provided with more 
befitting workout plans. Academicians should educate young students about 
PCOS, its risk factors, causes, prognosis, symptoms, diagnosis, and 
available treatment options. These efforts can make early detection of 
PCOS more likely. 


Women who have a BMI above 23 should be informed about its dangers 
and risks. They should be advised to lose weight in order to maintain a 
healthy body and mind. Women exhibiting signs of hyper-androgenism or 
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menstrual irregularities should be checked up to ensure the early detection 
of the syndrome. Timely identification of the disease combined with 
immediate therapy enhances the quality of life and helps to avoid health 
risks. This study can be used to explore the obstructions women face in 
lifestyle improvements such as regular physical exercise, weight control, 
following a healthy diet plan, and keeping their menstrual cycle in check. 
This survey can also be used to observe the effects of these improvements 
in a test group of women. 
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ABSTRACT 


Geometric modelling of several intricate and complex structures such as 
trees, mountains, clouds, ferns, geographic topography, and coastlines is 
challenging in computer graphics. Traditional splines such as trigonometric, 
polynomial, exponential, and rational fail to simulate this significant class 
of complex structures, which are highly irregular in nature. For this purpose, 
this research develops a novel cutting-edge method for synthesizing and 
modelling structures. The proposed technique; C! fractal interpolation 
function (FIF) builds an iterated function system (IFS) by integrating fractal 
calculus and rational cubic polynomial functions. Appropriate conditions 
on scaling and shape parameters are derived to help maintain the inherited 
shape qualities of the data. Experiments in numerous scientific domains, 
such as the pharmaceutical and chemical industries have been presented as 
an example, to confirm the usefulness of the suggested model. Moreover, 
the graphic results demonstrated that the developed monotone hybrid model 
(MHM) offers a heterogeneous method for gathering data with a monotone 
structure. 


Keywords: fractal interpolation, hybrid model, monotone interpolation, 
spline functions 


1. INTRODUCTION 


Interpolation is an indispensable tool for estimating the unknown 
function for data modelling. There are numerous interpolation techniques, 
which use various function families, including polynomial, exponential, 
rational, trigonometric, and spline. However, these traditional non-recursive 
approaches resulted in interpolants, which were differentiable several times, 
with the possible exception of a limited number of points and restricted 
model smooth structures. Many signals in the real world and in 
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experimental data are complex and their composition frequently 
demonstrates irregularities and hardly ever appears to be smooth. Thereby, 
interpolants that could picture non-smooth signals, as well as smooth 
structures, were highly required. The prime purpose of this research is to 
provide an innovative framework for the approximate representation of 
non-smooth data from natural (feathers, trees, leaves, clouds, feathers, 
flowers, glaciers, landscapes, and galaxies), engineering or scientific 
phenomena together with the collection of smooth data. A novel form of 
fractal interpolation function is created utilizing the rational spline 
containing a cubic polynomial in the numerator as well as the denominator. 


One of the pivotal characteristics of data generated from countless 
physical phenomena, scientific problems, and engineering applications is 
monotonicity where entities only have a meaning when their values are 
monotone. For instance, in the field of engineering, Tensile Strength, 
Ultimate Tensile Strength or Ultimate Strength is not an uncommon 
concept. It refers to the maximum level of stress that a structure or material 
can tolerate while being pulled or stretched before getting breaking. This 
phenomenon gave rise to monotone data. The applied forces, commonly 
called stress and the stretch, referred to as strain, are always in a monotonic 
relationship where one varies with the change in the value of the other [1]. 
Certain devices, such as; digital-to-analog (DAC used in audio video 
devices) and analog-to-digital (ADC used in music recording and digital 
signal processing), have monotonicity an essential characteristic used in 
feedback control loops [1]. In CAD and DAC, the digital output code 
always increases as the analog voltage input increase and vice versa. 
Monotonicity can be seen in Newton's law of cooling in which the rate of 
heat loss is directly proportional to the differences in temperature between 
the bodies [2]. 


Monotonicity also appears in other phenomena, such as the content of 
atrazine and nitrate in shallow ground waters, the link between the partial 
pressure of oxygen, and the percentage dissociation of hemoglobin 
percentage dissociation and oxygen partial pressure, erythrocyte 
sedimentation rate in patients with cancer and concentration of uric acid 
level in a patient affected from gout [3]. It is a fundamental function behind 
the renowned ‘Probability theory’ that uses the monotonically increasing 
function of a random variable in its cumulative distribution. 
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Natural phenomena can produce monotonic positive data as well. The 
relationship between a child's height and their intake of food and energy can 
be considered as an example. For instance, height is a positive quantity and 
always goes up. There are several other examples in everyday experiences 
where objects only create meaning when their values are monotonous. 
Examples are surfaces and the dose-response curve in biochemistry and 
pharmacology; consumption function that describes the relationship 
between a household's consumption and its available income in economics, 
statistical approximations of copulas and quasi-copulas, [4] and empirical 
option pricing models in finance [4]. 


Based on different approaches (such as classical spline, spline with 
knots, and fractal spline) various monotonicity-preserving models for 
monotone data have been discussed in the previous literature [1-18]. Spline 
functions have also been used for the solution of partial differential 
equations [19, 20]. Previously, spline functions have become the main tools 
for solving monotonicity problems. For instance, Crisp and Hussain [5] 
solved the issue of monotonicity through Bernstein-Bezier rational cubic 
functions. Spline functions were also used to provide multiple kinds of 
monotonicity methods for bivariate functions on triangles [6]. 


Using different types of spline functions, some other useful techniques 
and algorithms were supplied to address the problem of monotonicity [1, 3, 
7]. Typically, data that comes from experiments (for instance, information 
drawn from signals in the actual world, including seismic information, 
financial series, and bioelectric recordings are complex and immensely 
irregular. As a result, interpolating these data using traditional interpolation 
techniques became an inappropriate challenge. Barnsley [8] made an 
attempt by introducing a new interpolation method called Fractal 
Interpolation using a special type of iterated function system. It provided a 
powerful framework to deal with highly irregular data. However, 
consistently, when it is a prerequisite that the interpolant should reflect the 
inherited shape features (convexity, positivity, monotonicity) of the 
irregular data set, fractal interpolation functions failed to secure the shape 
properties of data. Later, Barnsley and Harrington [9] noticed that in case 
the problem is of shape preservation type, then the parameters of the IFS 
may be chosen suitably so that the corresponding FIF preserve the inherited 
properties of the data. This observation initiated a striking relationship 
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between the classical splines interpolation functions and the fractal 
interpolation functions. 


Recently, with the help of Barnsley and Harrington's results, various 
spline functions have been generalized in the form of fractal spline functions 
for instance [4, 10]. Fractal spline [11] though interpolating both regular 
and irregular data does not guarantee that the model satisfies the monotone 
trend lies in their data. Non-monotone visual models of the aforementioned 
experiments misguide us about the actual scenario (which might help 
scientists to recognize and evaluate various patterns and artifacts in their 
data) in these experiments. Hence, the visual model must maintain the 
geometric property of the experimental data in order to achieve precise and 
proficient results. 


The goal of the underlying research is to build a fractal geometry-based 
function to simulate hidden patterns and future behaviour in many 
disciplines. A monotone hybrid model (MHM) which is a combination of 
fractal interpolation function (FIF) and spline interpolation model (SIM) 
was presented for this purpose. The shape-preserving model is then created 
with the help of the proposed function to handle the non-linearity problem 
that the traditional interpolation model was unable to solve. Additionally, 
the proposed class surpasses its conventional non-recursive counterpart in 
approximating monotone functions with varied irregularities/fractality in 
their first derivatives (smooth to nowhere differentiable). 


The proposed hybrid model offered a variety of shape control 
parameters that can be sufficient and highly useful to achieve desired 
results. To meet the objective, a general piecewise rational cubic fractal 
interpolation function (GPRC FIF) [11] is used in the construction of MHM. 


The proposed model has the following useful features in comparison 
with the existing models: 


e The proposed scheme converts the hybrid model [11] to a monotone 
hybrid model. 


e The proposed hybrid model MHM is the generalization of the model 
presented in [7]. 


e Unlike [12], no derivative constraints are imposed in proposed model 
as MHM works for both data and data with derivatives proficiently. 
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e Once, the shape parameters are selected, the curve representation 
becomes unique in its description. 


e Existing methods have no degree of freedom [1, 3, 12, 13], whereas 
the proposed model gives two degrees of freedom. 


e The proposed model works for larger as well as smaller data, whereas 
in [14], constraints are developed on the subinterval lengths for data 
modelling. 


Four shape parameters are used in the proposed model, whereas in [3, 
15] one and two shape parameters are used respectively. The remaining 
sections are structured as follows. Section 2 focuses on a brief introduction 
to the basic methodology of the hybrid model developed in [11]. In 
Section 3, a new scheme is developed, which converts the hybrid model to 
a monotone hybrid model, for solving the problem of monotonicity. 
Whereas the proposed model is tested on different real data sets to verify 
their usefulness in Section 4. Comparative analysis has been presented in 
Section 5 and finally, Section 6 concludes with the significance of the 
proposed study. 


2. METHODOLOGY 


In this section, the hybrid model (GPRC FIF) [11] based on the fractal 
framework and spline function was briefly reviewed. The fractal framework 
was built on an iterative function system (IFS), which is defined below 
comprehensively. 


2.1. IFS Theory 


Let the real interval I[r;, m] be partitioned such that 7, < r; < = < %. 
Suppose the data set is given ((r, 9j) € I X E:i = 1,2, ..., n), where E is 
denoted as a suitable compact set containing all $j's. Take I; = [7,744], 
then the two mappings given below would be contraction homeomorphic so 
that €i: I Ij, 


with 
&gn)2rn, &(n)-2r,,i-212..n-1. (2.1) 
le;(e,) — &(e3)] < lile, — e2|, Veg,e; eL fosome 0<1; <1. 


Suppose, A = I X E and there is a mapping ci: A > E, which is 
continuous if 


2 — SIR 


Scientific Inquiry and Review 


Volume 7 Issue 1, 2023 


Arooj et al. 


a; (11,01) = Vis Oi (Taw Vn) = Visas T= 1,2, sag ll = 1. (2.2) 


lo;i(r,x) ^ o;i(ry)| € xilx - yl,re 1, Vx,y € E, forsome-1« kj < 
1. 


Take Tj:À > A, which is defined as: 


T(7,0) = (sir), a;r, $))v (r,$) € A. (2.3) 
Where 

&(r) = a;r + bi, o;(r, $) = kið + gi(T), 

with a, = ft, p, = te 


Here, gi(r) is any suitable continuous function. The set {A;T;:i € 
(1,23, =, ,n — 1)} is termed as an iterative function system. 


2.2.Hybrid Model 


The hybrid model [11] q is termed as GPRC FIF corresponding to the 
iterative function system as: 


p(e(r)) = kiel) + g(r), (2.4) 
- — Pir) _ Pi(R) 

with gi) = Oey = Y 

where 


P,(À) = ni(zi — ki$3)(1 — A)? + (Ni hidi — On — rnikids + ôi; — 
Ki0,)) 8 (1 — Â)? + (—ajhidisa + On — r1 )wikidn + pii — 
K;i$,)) ^ (1 — Â) + Wiig. — kid) P, 

Qi(u) =n - A) +i A- AY * pil A- R8) * eil?, â = 
L3 j$e€[nr. 


Mm-Ty 


Here, k; < 1 are the scaling factors, d; is the derivative value at the knot 
points and ni, Ôi, pj ,@; are the non-negative shape parameters. 


3. A MONOTONE HYBRID MODEL 


This section sets up a monotone hybrid model (MHM) to solve the 
monotonicity problem present in the hybrid model as defined in [11]. In 
setting up the MHM, certain conditions were imposed to shape the 
parameters and scaling factors, to acquire the desired model. 
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Let, the monotone data set {(17;,0;) € I X E:i = 12, ..., n} be distributed 
such that r € r; €: € r, with Ùi 2 Ùi, i= 1,22, ...,n — 1. 


Or alternatively, 
238—220, (-12,..n-1. (3.1) 
i 
For the function that increases monotonically, the necessary constraints 
were applied on derivative as: 
ði = 0,i = 1,2,...,n. (3.2) 
For the data set increasing monotonically, there are two possibilities. 


Case 1: If č; = 0 and k; = 0 particularly for some monotone and collinear 
data set, the function defined in Eq. (2.4) is constant in each subinterval. 
Therefore, 


0; = Oind = 0,i = 1,2, sse y Tos 
So that 
o(r) =%;, Vvrelryrl]i-212..,n- 1. 


This demonstrates the process by which a hybrid model automatically 
transforms into a monotonous hybrid model. 


Case 2: If č; > 0and č; — kičn > 0e x < Tuc , Where &, = 
n "1 


mE “1 for this case, we have the following theorem. 


l 


Theorem 3.1: Suppose {(r;, i) : i = 1,2, :-:, n) be a monotone data set 
such that $; > 0. Assuming n; > 0, w; > 0 and ð; constitute the derivative 
at the knot points ?;. Then the corresponding hybrid function GPRC FIF 
converts to monotone hybrid model (preserve monotonicity) in each 
interval if the scaling factors and shape parameters satisfy the following 
condition: 


(i) The scaling factors k;, i = 1,2, ..., n, are picked as: 
H * * * 
0 € k; < miníai, ki; K5j, K3;] 


(ii) The shape parameters 6;, pj, i = 1,2, ...,n, are chosen as: 
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K; K; 
wi(dix1-G4n) (meee J 
AM. .QEOAEC TT y SS aa . 


> ——_—, pi > max i 
(roe) Pi | Gio) 8ieme ai= zta) 


where 
* ai Oi41 * aiði x _ Visi-Vi 
1i ðn > 2i 04 > 3i 94—94 


Proof: Take the monotone data set into consideration without losing 
generality and 


e(e;(r)) = kjQ(r) + gi(r) 


It is clear from single variable calculus that the function, «(ei (r)) is 
monotone increasing (decreasing) if and only if pY > 0 (@Y < 0). 


Here, «p is FIF, so derivative qO is also FIF and satisfies functional 
equation. After careful calculation, we arrived at the following equation by 
differentiating the above equation with respect to r. 


(^ ene 
e(uG)) 2 "2 g (r) (3.3) 


The hybrid function defined in (4) becomes monotone hybrid model if 
qO > 0 and only if g; = 0 for all r e [r,rj,4]. 


Taking that x; > 0 v i € {1,2,3,---, n], itis evident that pY 2 0,vr € 
[ri ri if 


a Bean y 
gi? (r) = X46; =E z0,vue [04]. (3.4) 
[ai(& )] 
where 
Gi = wf (Oii — En), 


€5; = 6g + Za (01; — kin) — ni(di — 791); 

C31 = (63; — ĉii) — (30a; — &ipi (i — KiEn)) + (6,0; (ði — =, On) F 
nip: (ð: — 1:84). 

Cai = (Csi — ĉesi) — (3Niwi — Sipi (Ei — Kin) ) + (Siwi (Oi41 — PL + 
nipi (ði — 7-91); 
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Cs; = Ĉesi + 2A (pii — ibn) — wiit — Fn); 
6, =7.2(9, -“a 

Cei = Ni ( i a; 1)- 


The denominator in Eq. 3.4 being a squared quantity is always positive. 
Furthermore, if 9; 20 and w;=0, the sufficient conditions for 
monotonicity in each subinterval [r;, r;,.4] are: 


x x Ai Oi+1 

Gi 20, j-1,2,7,6 ĉi 20S ki S 7, (3.5) 
n 
6,20 ek S E (3.6) 
1 

Also from Eq. (2.1), it implies that o € C*[1r,, m], 
whenever 
—aj < K; < a; 
and 

& —k£, > OK < i (3.7) 

Vn-V1 
However, ¢2; >0 & 6; > 

n(4—121) 

(izkin) ` 2) 

a x o(05 71x) 
Similarly, Csi > 0S pi > eae (3.9) 
Finally, both €3; > 0 and &4; > 0 iff 

Sici -Bi9i(0141 7 an) 

> —l , 3.10 
"i 8, G7) mi (477101) Pu 
provided 

n(2i-221) 

A eI" 


Hence, the results of Eq. 3.5-3.10, proved the theorem. 


Corollary 1: If x; =0,i=1,2,..,n—1 in above Theorem 3.1, the 
sufficient condition convert to 
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nii Wiði+ı  (3NiWi—SiWi9i41 
Bi ü Pi > max f Ši ( oisi *nii )} 

This is sufficient condition of monotonicity for classical GPRC. 
4. RESULTS AND DISCUSSION 


This section presents two experiments, which demonstrated the validity 
and adaptability of the proposed MHM. In these experiments, the 
derivatives values were calculated through geometric mean method. The 
monotone data set taken in Table 1 presents the experimental result 
concerning the watering the great northern beans. Distilled water with pH 
of 8.5 and Potassium Hydroxide was combined in the chemical solution, 
which was later used to watering the beans. After 40 days, the beans plant 
was taken out from vat to weight for the effect of chemical solution. It can 
be easily noticed that resulting data is monotone. The r-values represented 
the days and the v-values showed the height of beans. 


Table 1. Height of Great Northern Beans 


i 1 2 3 4 5 6 7 8 
ri 1 2 12 18 24 30 36 40 
ù 0 0 042 2.08 3.43 3.78 4.12 4.37 


The data demonstrated in Table 1 is monotone; therefore, the resulting 
curve must be monotone. However, from Figure 1, which is constructed 
through hybrid model [11] it can be observed that the arbitrary value (7;=7, 
ôi=4, p;—8.2, «0;-99, K;=0.9) does not preserve the inherited monotonicity 
of the curve. On the other hand, Figure 2 displays multiple curves, which 
are produced through the proposed MHM. In Figure 2, it is found that the 
curve at different iterations preserve the essential characteristic of the data, 
which is monotonicity, proficiently. Figures 3 and Figure 4 highlighted the 
effect of diverse values of A; and œ; ,rescpectively. The data presented in 
these figures showed that with different choices of parameter values through 
MHM, monotone curve of any particular requirement may be easily gained. 


The monotone data set in Table 2 reports the varying degree of 
saturation of hemoglobin with partial pressure of oxygen. It was observed 
that saturation of hemoglobin was directly proportional to the partial 
pressure of oxygen. The r-values represented the partial pressure of oxygen 
and z-values indicated the percentage of saturation. 
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Table 2. The varying ability of hemoglobin 


i 1 2 3 4 5 
" 0 2 8 10 18 
$; 0 70 91 91 110 


The data demonstrated in Table 2 shows the monotone relationship, 
therefore, the curve must exhibit the same characteristic. Figure 5 is 
generated through the random values to the shape parameters and scaling 
factors (1;22.9, 6;=3, pi=7, wi=15, k;=0.7) in the description of hybrid 
model as defined in Eq. (2.4). Figure 5 depicts non-monotone hemoglobin 
dissociation curve (HDC) which does not make any sense, as the behaviour 
of the curve should retain monotonicity. This flaw is removed in Figure 6 
by using the proposed method presented in Section 3. It can easily be 
observed that Figure 6 conserve monotonicity at each level of iterations. 
Figures 6 and 7 showed the profile of HDC with diverse values of n; and cj. 
It is worth noticing that although various parameters may lead to different 
results but all conserve monotonicity. 


Height(cm) 
Height(cm) 


9 Data 

77 (st iteration 
2nd iteration 
piang 3rd iteration 


“0 5 10 5 0 2 3350 0 5 10 15 2 25 30 35 40 


Figure 1. Hybrid model [1] Figure 2. MHM using 7; = 0.5, 
Qj = 70 
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Figure 5. Hybrid model[1] 
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Figure 7. Effect of 7; using 


Qj = 3.8 
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Figure 6. MHM using 7]; = 
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4.1. Comparative Analysis 


The proposed method guarantees to produce monotone model as 
compared to the hybrid model [11], which does not ensure to exhibit 
monotone model when data is monotone. All the existing models [3, 12, 
13], failed to provide even a single degree of freedom, so MHM was a 
noteworthy addition as it offered two degree of freedom to modify and 
enhance the model. Unlike, [14] and [12], it is easy in the applicability, as 
it does not constrain the interval length or impose any conditions on 
derivatives, respectively. The scheme developed in [7] only works for 
regular data but the proposed approach, which is generalization of [7] is 
applicable for both regular and irregular datasets. 


4.2. Conclusion 


General Piecewise Rational Cubic Fractal Interpolation Function 
(GPRCFIF) was developed to provide a deeper understanding of the data, 
which displays rough, irregular, and fragmented configuration. Generally, 
it originates from complex functions or experiments. The GPRCFIF offers 
the family of four parameters and one scaling factor to modify, upgrade, 
and change the visual images. One of the most prominent characteristics of 
the GPRCFIF is its capacity to visualize data with smooth curves by setting 
all scaling factors to zero. The current study, in particular addressed the 
problem of monotonicity through the novel proposed model MHM. The 
experimental results illustrated that the performance of the proposed MHM 
was very efficient and heterogeneous as it contains the ability (due to 
recursive nature) to produce multiple monotone curves for the same 
monotone data according to the user's desire. 
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ABSTRACT 


A powerful tool for dealing with ambiguity, attributive values, fuzziness, 
and inconsistency is the cubic intuitionistic fuzzy soft set. In this paper, 
multi-criteria group decision-making (MCGDM) technique called extended 
TOPSIS based on distance and similarity measures is presented. 
Furthermore, a real-world problem is resolved by applying CIFS-set to 
support the theoretical visualization. Arguably, the proposed technique 
works well and has practical uses. 


Keywords: cubic set, intuitionistic fuzzy set, intuitionistic fuzzy soft set, 
multi-criteria group decision-making (MCGDM), TOPSIS technique 


1. INTRODUCTION 


Taking the best decision requires taking into account a range of factors 
due to increasing complexities in business, engineering, scientific, and 
technological environments. The most effective way to select the ideal 
alternative among all possible options is decision-making. Since 
generalized variants are frequently used in decision-making, nearly every 
other problem has a considerable range of requirements. Such requirements 
frequently come into conflict with each other and no solution may ever be 
able to fully satisfy all of them. To address such challenges, decision- 
makers need to overcome the MCGDM problem. For this purpose, many 
mathematical theories have been expounded such as the fuzzy set theory, 
intuitionistic set theory, and interval-valued intuitionistic set theory. In 
1965, fuzzy set theory was introduced by Zadeh [1]. This theory deals with 
membership value over closed interval [0,1] and overcomes the vague and 
ambiguous environment. However, due to the increased complexities in the 
environment, a decision-maker faces difficulties to express their idea in the 
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form of a single membership value. To deal with such a situation, the 
intuitionistic fuzzy set (IFS) [2] theory was introduced. This theory not only 
deals with the membership value but also with the non-membership value 
and hence gives a more precise result which helps in decision-making. In 
this context, many researchers have presented different theories based on 
distance and similarities measures [3—5]. 


The technique known as the order of preference by similarity to ideal 
solution (TOPSIS) was introduced by Lai et al. [6] and remains a widely 
known method. The purpose of this method is to identify the longest path 
from the negative ideal alternative (NIA) as well as the shortest path from 
the positive ideal alternative (PIA). 


Since this technique has been introduced, various researchers have used 
the TOPSIS method to resolve decision-making problems in an ambiguous 
environment. Interval-valued intuitionistic fuzzy set and its fundamental 
properties was proposed by [7]. Jahanshahloo et al. [8] presented the 
TOPSIS method for DM under a fuzzy environment, Chu ef al. [9] 
suggested the TOPSIS method for robot selection, and Zulqarnain et al. [10] 
generalized the fuzzy TOPSIS for MCGDM. Saglain et al. [11] discussed 
the use of the TOPSIS method for the selection of a smartphone. Shen et al. 
[12] extended the application of this technique under intuitionistic fuzzy 
environment, Chang et al. [13] discussed the distance approaches using the 
TOPSIS method, Li et al. [14] suggested the Pythagorean fuzzy TOPSIS 
based on similarity measures, and Gupta et al. [15] discussed the extended 
TOPSIS under interval-valued intuitionistic fuzzy environment. 


All these theories involve uncertainty to some extent as they do not deal 
with the membership and non-membership values simultaneously along 
with fuzziness. To deal with such kind of situations, Garg [16] proposed the 
TOPSIS method for MCGDM under a cubic environment which deals with 
both the intuitionistic value and fuzziness, simultaneously. Pramanik et al. 
[17] defined the TOPSIS method for neutrosophic cubic information, Jun et 
al. [18] suggested the cubic fuzzy set (CFS), while Garg [19] discussed the 
cubic intuitionistic fuzzy set (CIFS) and its fundamental properties. The soft 
set theory with its fundamentals was proposed by Maji et al. [20]. The 
theory of soft set was merged with cubic set and is known as cubic 
intuitionistic fuzzy soft set (CIFSS). It was proposed by Saglain et al. [21]. 


oe 
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Cubic intuitionistic fuzzy set (CIFS) is considered highly effective for 
decision-making problems as it involves interval-valued intuitionistic fuzzy 
number (IVIFS) and intuitionistic fuzzy number (IFS) over the interval. 
Still, this theory does not deal with multi-attributive values. Considering the 
soft set, this paper attempts to define the multi-attributive decision-making 
problem under the cubic intuitionistic fuzzy soft environment, where each 
element is defined by the CIF-number (CIFN). A methodology which 
utilizes the extended TOPSIS method is also proposed. Furthermore, some 
distance and similarity measure formulas have been defined which evaluate 
the positive ideal alternative (PIA) and negative ideal alternative (NIA). 


2. PRELIMINARY SECTION 


In this section, some basic concepts of intuitionistic fuzzy set (IF-set), 
interval-valued intuitionistic fuzzy set (IVIF-set), soft set, intuitionistic 
fuzzy soft set (IFS-set), CF-set, CIF-set, and CIFS-set are defined. 


Definition 2.1 [2]. An intuitionistic fuzzy set (IF-set) K over a crisp set U is 
defined by a function 1j; represented by the mapping 


w = U > [0,1] x [0,1], 


where qy is the membership function which deals with truthiness and 
falsity and U be a universal discourse. Thus, an intuitionistic fuzzy set (IF- 
set) over set U can be represented as 


k= {(3, (t« C, fs C2) | eU] such that 0 € tg(z) + fa(z) € 1; VzeU. 


Definition 2.2 [20]. Let U be a universal discourse. An interval-valued 
intuitionistic fuzzy set (IVIF-set) K over a crisp set U is defined by a 
function 1j; represented by the mapping 


u = U > P([0,1] x [0,1]), 


where y; is the membership function which represents the degree of 
belongingness and the non-belongingness value. Thus, an interval-valued 
intuitionistic fuzzy set (IVIF-set) over set U can be represented as 


k = (s, CI C), 6 C2 7, [f C, f C) DAEL} such that 0 
X tex) + fa(3) € 1; VaeU. 


Definition 2.3 [19]. Consider U to be a universal set and € = 
(€ Cz, C3, ...., Cn} be the set of attributive values. Then, a soft set 4 over a 
crisp set U is defined by a function øg represented by the mapping 
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Qgy C P(U), 


where qs is an approximate function and its value qs(C;) can be 
represented as 


§ = f(c PED) ici € Gos) € PO}. 


Definition 2.4 [17]. Let U be the set of universal discourse and 3 be a cubic 
fuzzy set (CF-set) over U and defined as 


3 = UC, WG), v) ZEUS, 

where ug(3) = [Ug(3)~,Ug(v)*] and v represents the interval-valued 
fuzzy set and FN in 2€U, respectively. 
This pair is denoted as 3 = (Ug, v), which is known as a cubic fuzzy set. 
Definition 2.5 [18]. The term cubic intuitionistic fuzzy set (CIF-set) 3 
defined over U is defined as 
3 = {(C, Ue), v)IaeU], 


where "uk (3) = UC, (It (a) be C0 * T Df C2), f C2 D, and v = 
[tg Cz), fa (2)] represent the interval-valued IFS-set and intuitionistic set (IS) 
respectively such that 0 € tg(3) +fg(3) < 1. Also, 0 € t) fa() < 1. 


This pair is denoted as 7 = (1j, v) and termed as a cubic intuitionistic fuzzy 


set. 

Definition 2.6 [18]. The term cubic intuitionistic fuzzy soft set (CIFS-set) 
3 defined over U is given as 

3 = {cy wGi), vils;eU), 


where (ai) = (ci (BYS tt] Gri) RBD I and vi = 
[t Cai), fa Cui] represents the interval-valued IFS-set and intuitionistic set 
(IS) respectively such that 0 < tg(3:) + fai) € 1.Also,0 € 
be (3a), fa (3i) S 1. 
This pair is denoted as 3 = (uj, Vi) and termed as cubic intuitionistic fuzzy 
soft set. 


Definition2.7 [18]. Consider 
B= (eG), 55 327], eee 
([ts i), fs i) ]) 
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(ty Cg); fy Cu) 


Then, for q = 1, the distance measures are defined as 


y- [nece ected — M 


e Distance measures 
dq(B,¥) = (Xh tlte G) — by (Si) 14 + Ies (32* — by CR) * I? + 
fai” — fyi I? = BBD — fy(30* I? + te — ty (301? + 


la(s) -KEDD 

e Normalized distance measures 

dq'(B,¥) = (LEP sflts(i)7 — ty(&) |" + [bs G@i)* — &y(30* ]* + 
IBD — fyi I + Ifa GD" — fy 30" I? + Its (a) — ty (321? + 
la(s) —f&(019))* . 

where n represents the mean of attributive values. 


3. EXTENDED TOPSIS TECHNIQUE 


In this section, based on the suggested distance measure, a TOPSIS 
(technique for order of preference by similarity to ideal solution) technique 
for tackling MAGDM problems in the form of CIFS sets is proposed. 


3.1 Description of the Problem 


A recently opened restaurant is hiring a new head chef to take charge of 
the kitchen. For hiring a head chief, they have published the advertisement 
in the newspaper and different applicants have applied in response. Assume 
that there is a set of m applicants (alternative) 3; = (41, 32,33, =- &m) 
chosen for the interview. The restaurant has gathered decision-makers 
(T1, T?, T?,T^) and assigned them the duty to identify the best head chef for 
the restaurant. The selection committee has decided to assess the applicants 
% = {41,32 52 +» 3m} based on n different criteria d = (di, do, d3, ..., dy}. 
All applicants participated in group talk for this purpose and a panel then 
developed results for each applicant in the form of CIFS-set criteria X;; — 


(wij, Vij), where 


mij (Bi) = (EuB bi BD L ij SD fy GS *]) and vi; = 
[bij (i), fij (4))] represent the interval-valued IFS-set and intuitionistic set 
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(IS), respectively. Here, the components [t; jGD bij (3;)*] and t; ne 
show that satisfaction degree corresponds to a given alternative H; which 
satisfies the criterion Cj, whereas the components 
fi 30 ag BD] and fi; (3i) represent the degree of dissatisfaction with 
the given alternative 3; regarding the criterion dj. Thus, overall results are 
recorded in the form of a CIFS-set environment and defined in matrix form 
which is represented as D = (X) mxn- 


3.1.1. Computing CIFS-PIA and CIFS-NIA. The CIFS-positive ideal 
alternative (CIFS-PIA) and CIFS-negative ideal alternative (CIFS-NIA) 
regarding the alternative z;(i = 1,2,...,m) may be selected as 1 and 0, 
respectively. It shows the ranking values of each alternative. Thus, the 
ranking values of CIFS-PIA and CIFS-NIA can be represented as 


* = (([1,1], [0,0], (0,11)4«4 and X^ = (([0,0], [1,1], (1,0). 
So, X* and X~ complement each other. 


If the decision-maker wants to define the other reference points as 


= = (dt "2d DT 25 £P n (p*,q +), (3.1) 
and 
X^ = (E LGA) erar) (3.2) 
where 
tht = max(ti (s) ), tj = max(tij(4i)"), fj = min(fij (i), p 
= min(f;; G)*), 


pj — min (ti (30), qj = max (fij(3i)) vi 

and tj- = min(ti (s) ) Qj o- min(ti; (s))*) dj 
max(fu (s) ), ff = max(fij(3i)"), pj = max (tua) q 
minj (WED) vi 

3.1.2. Computing Distance Measures Between Alternatives. 
Consider the criteria in terms of the weight vector w; = (w1, W2, ..., Wa) 
along with CIFS-PIA(X*) and CIFS-NIA(X ). We evaluate the weighted 
distances between the alternatives 3; and its X* as well as X^ defined as 
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dg (ai, X*) = | yum w;{|t;* = tug] + [trt = ij (3i)* |" T 


ICONE fe + |f BD+- m + |b) — pt |" +|- 
1 
/ 
EDY p (3.3) 
and 
da (3,3 )— [: jam wil|ti; (i) = vali + |i; C2* E "mu t | fj lk 
1 
m z _ / 
f(a)" |" +f —f@i)* |" + |p; - a" + |f; 2 - 4j Ml - 
(3.4) 
whereq 2 1 € R”. 


The relative closeness coefficient of alternative z;(i = 1,2, ..., m) with 
respect to CIFS-PIA (3*) (which is based on weighted distances) is defined 
as 


o dq (5X ) qo. yt 

LU dG GX) +dq GXH)” dq (3i¥") *0. (33) 
Further, 

0 < dg (3i X) < dg (3X) + dg, X*) and 0 < C; < 1. 


3.1.3. Group Decision-making TOPSIS Approach. Using the above 
analysis, an approach is presented to solve the group decision-making 
problem under the CIFS-set environment. For this purpose, consider there 
are 'T' decision-makers {T+,T?,...,7*} which evaluate the set of ‘m 
alternatives z;(i = 1,2,...,m) under the set of ‘n’ different criteria d = 
(di, dz, d3, ..., dn}. 


These decision-makers give the results in the form of CIFS-set (Xi Pu 
(Cay (vis) ), 
where (ij; GD)“ 
= CRODH yB), [dijo )*, Gu G0*)*]) and (vij)* 
= (ti) i GOD", where K = 1,2, K. 


Further, suppose that the criteria of the weight vector (w;)* = 
((w1)¥, (w3)F, ..., (w4)*)? such that each (w;)* > 0 and X? 4(wj)* = 1. 
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Moreover, to overcome the diverse judgments by different decision-makers, 
their decision is prioritized according to weight vector y= 
(Y, Y2, - Yi), such that y, > 0 and XE 4 y, = 1. 


The stepwise algorithm of TOPSIS is presented below. 


Step 1: Arrange the ranking values of the alternative given by each 
decision-maker in the form of a matrix 


d, di oo. d, 
- - XQ Gu a. x0 

k) 
T EE Xa X287 Zr Kea” 
(m : : ^ : 
La Je uu See 


Step 2: For each decision-maker T* = {T1,T?, ..., TK), compute CIFS- 
PIA and CIFS-NIA corresponding to the applicants z;(i = 1,2, ..., m) using 
the Eq. (3.4) and (3.5) respectively and defined as 


Ot) = (Ee e» LG GP) (eps ap) G9 
and 
OF = (IG GY LLP, PY) (ep Gi), = B.D 
where 
(ti*)* = max{ (ti (i) 269], 07^ = max((5u(30*)09], 7)" 
= mint Gug 5)" = min(Qi(a0*)99], 
(pi)* = min((30)9], (qj)* = max(; (009) Vi 
and 
(t )* = mint (bij) )09j, (t) = min(Q(30*)9], Gy 
= max(ij (a) 99], jp max{ (fij(3i)*)}, 
(pj = max{(t5 (30) ha 
= min; TOM vi. 
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Step 3: For each decision-maker, compute the separation measures 
between the alternatives z; from its CIFS-PIA(X*) and CIFS-NIA(X ) 
which are denoted by dj ((z;)*, (X*)') and d; ((3:)*, (X )*), respectively. 


Step 4: For each decision-maker, the relative closeness coefficient is 
determined as 


dg (GO*, 0:9) ag (GDEK) k = 1,2,...,K (3.8) 


where d4 (()*, (S*)*) # 0. 


c; = 


Step 5: Rank the alternatives based on the descending values of Cjs 


Group decision making TOPSIS approach 


O 2 compute CiFS-PIA and CIFS-NIA 


0 3 Compute separation measure 


0 4 Determine relative closeness 
coefficient 


05 Rank the alternatives 


Figure 1. TOPSIS Algorithm 


Example: To demonstrate the above-mentioned approach depicted in 
Figure 1, an example is discussed below. 


4. RESULTS 
4.1 Case Study 


A recently opened restaurant is hiring a new head chef to take charge of 
the kitchen. For hiring a head chef, they have published the advertisement 
in the newspaper and different applicants have applied in response. A total 
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of four individuals 3;; i = 1,2,3,4 have been chosen for the interview. The 
restaurant has gathered decision-makers (T1, T?,7?, T^) and given them the 
responsibility to identify the best head chef for the restaurant. The selection 
committee has decided to assess the applicants z;; i = 1,2,3,4 based on four 
criteria d = (d, d5, d3, d4} defined as 
dı: cooking experience d,: time management, dz: presentable, d4: team handling. 
For the assessment, they conducted group discussion with all the applicants 
and the results are formulated by a panel in the form of an IVIFS set. From 
lots of applicants appearing for group discussion, only four applicants are 
shortlisted for the interview. At this stage, the results are recorded in the 
form of an IFS set. Then, the following steps of the proposed approach are 
executed to find the best head chef for the kitchen. 
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4.1.1. Step 1: Ranking Values of Alternatives Corresponding to Each Decision-maker 


Table 1. Ranking Values of Alternatives Corresponding to T® Decision-maker 


Dacisi Applicants 
mn a a : 5 
eighs 
: ( [0.30,0.50], [0. aue ( [0.25,0.35], [0. Aona ( [0.32,0.40], [0. 2041) (s. 40,0.44], [0.50,0.52]), 
i , (0.35,0.29) (0.19,0.27) (0.35,0.49) (0.55,0.35) 
(€ [0.15,0.30], [0.35,0. 4 (a [0.15,0.18}, [0.20,0. 30h) (Paora 0550620 (© 55,0.59], [0.18,0.33]), 
x ,(0.17,0.40) (0.29,0.63) (0.40,0.40) (0.16,0.20) 
TO 3 (© 39,0.48], [0.10,0. 18h) (€ [0.23,0.32], [0.47,0. Sop.) eens) [020025] (s. 37,0.42], [0.40,0.50]), 
a (0.20,0.15) (0.30,0.60) (0.29,0.19) (0.30,0.70) 
(© 420.55], Ones Ga EU ( [0.16,0.32], [0. iam (© 20,0.55], [0.38,0.49]), 
n (0.15,0.30) (0.42,0.20) (0.17,0.34) (0.80,0.17) 
weights 0.40 0.20 0.25 0.15 
Table 2. Ranking Values of Alternatives Corresponding to T? Decision-maker 
2 di d, d; d, 
Weighs 
; ( [0.35,0.45], [0. sansa) ( [0.40,0.60], [0. 20 (© .16,0.35], 022,0320) ( [0.11,0.20] [0. Ce 
1 (0.60,0.27) ,(0.35,0.19) (0.32,0.20) (0.40,0.17) 
, (TS 150.30], [0.35,0. 49 (a [0.60,0.70], [0.19,0. 28h) (© .30,0.45], [0.21,0. 38h) Gas [o ozo 
z , (0.20,0.70) (0.60.0.10) (0.50,0.30) (0.40,0.30) 
TO : (a [0.19,0.39], [0.10,0. m (€ [0.20,0.35], [0.17,0. E (© .10,0.18], [0.12,0. 35h) at pns 
3 (0.11,0.20) (0.10,0.40) (0.20,0.35) (0.20,0.70) 
(pisos [0. [88A ( [0.15,0.40], [0. ET (e470) [032,0420 (esna anna) 
š (0.30,0.42) (0.30,0.63) (0.15,0.29) (0.66,0.18) 
weights 0.35 0.30 0.15 0.20 
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Table 3. Ranking Values of Alternatives Corresponding to T? Decision-maker 


a BAM d, d, d d, 
eighs 
, (© .30,0.50], Do Mon ( [0.25,0.35], [0. D ( [0.32,0.40], [038,0400 ( [0.40,0.44], [o 50,0520) 
: ,(0.35,0.29) (0.19,0.27) (0.35,0.49) (0.55,0.35) 
(© .15,0.30], on ( [0.15,0.18] [0. 2003911) ( [0.30,0.40], ssp] ( [0.55,0.59], [018033 
s2 , (0.17,0.40) (0.29,0.63) (0.40,0.40) (0.16,0.20) 
TO ü 39,0.48] [0. T ( [0.23,0.32] [0. supo (© 40,0.45], S (© 37,0.42], ea) 
5s (0.20,0.15) (0.30,0.60) (0.29,0.19) (0.30,0.70) 
Ga [0. Oen (© 44,0.50], s ( [0.16,0.32], 018022) (© .20,0.55], senem 
a (0.15,0.30) , (0.42,0.20) (0.17,0.34) (0.80,0.17) 
weights 0.40 0.20 0.25 0.15 
Table 4. Ranking Values of Alternatives Corresponding to T® Decision-maker 
Decision Applicants d d d d 
maker & Weighs t : 3 * 
(ito .20,0.35], ee) (s. 36,0.40], pps ( [0.30,0.55], [0. SAY ([0.11,0.20], [0.16,0.29]), 
ae (0.30,0.45) (0.45,0.20) (0.10,0.30) ( (0.35,0.19) ) 
à (te 180.30], [0.25,0. sena (qw. .18,0.30], [0.19,0. 3) (16150351 0260608) (1040020, [0.10,0.501) 
2 (0.20,0.70) (0.40,0.20) (0.25,0.40) (0.50,0.20) 
r6 (0 .10,0.29], [0.15,0. m (w .20,0.30], [0.11,0. 0) (Duos n (Dann sn 
ia (0.11,0.30) (0.30,0.50) (0.30,0.40) (0.25,0.60) 
(te 25,0.32], 400451. (‘eo .16,0.20], 02505405 ( [0.45,0.60], [0. aU ( [0.21,0.39], [0. 200498 
= (0.20,0.32) (0.10,0.63) (0.29,0.63) (0.66,0.28) 
weights 0.45 0.19 0.25 0.11 
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4.1.2. Step 2: Computing CIFS-PIA and CIFS-NIA for Each 
Decision-maker. Computing the CIFS-PIA and CIFS-NIA for each 
decision-maker is depicted below in Table 5. 


Using the Eq. (3.1) and (3.2), CIFS-PIA and CIFS-NIA are computed 
respectively as follows: 


(X* 9* = (Gr, GL [05 GP) (075 (729) 


and 


%7)" = ([c 5 e YL 105 0; F Iena) 
where 
(ty = max{ (ti (3) )}, (ay = max((5u(30*)09], FF)" 
= mint (fiji) }, FY" = mint (fiji) }, 
(pj) = min(( 009], (qj)* = max(G; 0009) vi 
and 


(t = min(Qu) J}. qr )* = mint (bij (%)")}, 7)" 
= max{ Gy (a) 2*9], A = max(u (3)*)}, 
(k) 
(Y = max f(t) pa 
(k) 
= minj (fuo) | Vi. 


Table 5. Positive and Negative Ideal for Each Decision-maker 


PIA 


NIA 
J [0.42,0.50], [0.44,0.50], [0.40,0.45], ((0.55,0.59], 
X. [fpo100.18]' (10.10,0.25]' (¥0.18,0.22]”” [0.18,0.33]’” 
T (0.17,0.40) (0.19,0.63) (0.17,0.49) (0.16,0.70) 
È 
[0.15,0.30], [0.15,0.18], [0.16,0.321, ((0.20,0.42}, 
(to350.45]" | | ‘fo.47,0.50)”’) [*[038049]" [0.50,0.52]"' 
1 (0.35,0.15) (0.42,0.20) (0.40,0.19) (0.80,0.17) 
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PIA 
DM 
NIA 
xt Goa ([0-60,0.70],, 047,050], (0-40,050], 
2 C 10.10,0.29]' [0.15,0.20] (10.12,0.15] /10.10,0.20]" 
m (0.11,0.70) (0.10,0.63) (0.15,0.35) (0.20,0.70) 
Te 
l0-10,0.29],, ({0-15,0.35],._ (0190-185, (0-11,020],, 
x- [0.48,0.54] [0.19,0.36] [0.32,0.42] [0.40,0.44] 
2 (0.60,0.20) (0.60,0.10) (0.50,0.20) (0.66,0.17) 
xt [0-30,0.50),, (0-30,0.38],_ TAE Rd (032,045],. 
3 [0.10,0.30] [0.20,0.30] 0.30,0.30] [0.10,0.20] 
T" (0.20,0.35) (0.11,0.40) (0.11,0.70) (0.30,0.40) 
TG 
[0100.35], (0-100-15],, (neuem ([0-10,0.29],, 
x- [0.40,0.45] [0.40,0.45] 0.50,0.50] [0.50,0.50] 
3 (0.45,0.10) (0.60,0.10) (0.25,0.32) (0.45,0.20) 
xt (0250351, in 030,055], [040,039], 
4 [0.15,0.35] [0.11,0.32]" (10.20,0.25] [0.10,0.29] 
(0.11,0.70) (0.10,0.63) (0.10,0.63) (0.25,0.60) 


(4) 


[0.10,0.29] [0.16,0.20] [0.10,0.29] [0.11,0.20] 
" o oa (aon ect quest d 
4 (0.30,0.30) (0.45,0.20) (0.30,0.30) (0.66,0.19) 


4.1.3. Step 3: Evaluation of Separation Measures. Without the loss 
of generality, assuming q=2 separation measures are calculated between the 
applicants from their CIFS-PIA and CIFS-NIA corresponding to each 
decision-maker and denoted by dz ((3i)*, (X*)*) and dj ((z)*, (X^ )*) 
respectively, as depicted in Table 6. 


Using Eqs. (3.6) and (3.7) 
dg (a, X^) = l: j=l w;{|t7* = tul t odd a tu)" [" 
BD- pi + |fij(ai)* - gs + |b) — pil" +|} - 
1 
/ 
A i 
and 
da 3i X ) = F j21 wil|ti;Cu2- x: a | + |i; C2* = Bs + | fj ian 
NU. U- 419 - q _yay)'/a 
HEDT + | f° -fyGOo*| 9|»; ^ ty Gl +I- 97 | j| 
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Solving these separation measures corresponding to each decision- 
maker, we get 


Table 6. Evaluation of Separation Measures 


TH TO TG TH 
Applicants [09 — q03 m e w Co 109 w 
k 0.3553 0.3102 0.4961 0.3927 0.2703 0.2501 0.3651 0.2461 
22 0.3927 0.4116 0.4087 0.4672 0.2180 0.2151 0.2583 0.2076 
33 0.3784 0.3478 0.2645 0.4783 0.2158 0.2426 0.2764 0.2123 
k 0.4067 0.3600 0.3538 0.3330 0.2058 0.2037 0.1947 0.1582 


4.1.4. Step 4: Computing the Relative Closeness Coefficient. Using 
Eq. (3.8), the relative closeness coefficient and the ranking order of the 
applicants corresponding to each decision-maker is computed as 


dol YD oL 
i dY (GF, EHE) + a! (GED)? k = 1,2,...,K 


where d4 ((s:)*, (X*)*) # 0. 
The outcome is manifested in Table 7. 


Table 7. Relative Closeness 


TD T TO TH 
Applicants c, "2 c, "o 
Z1 0.4661 0.4418 0.4806 0.4027 
Z2 0.5117 0.5334 0.4967 0.4456 
Z3 0.4789 0.6439 0.5292 0.4344 
Z4 0.4695 0.4849 0.4974 0.4483 


4.1.5. Step 5: Ranking of Alternatives. Using the relative closeness 
coefficient, the ranking of alternatives corresponding to each decision- 
maker is manifested in Table 8. 


Table 8. Ranking of Alternatives 


TO TO TG T 
Applicants 
cO Ranking c, Ranking cO Ranking cO Ranking 
Za 0.4661 4 0.4418 4 0.4806 4 0.4027 4 
32 0.5117 1 0.5334 2 0.4967 3 0.4456 2 
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TO TO TG T 
Applicants 
pp c, Ranking c, Ranking c, Ranking c, Ranking 
33 0.4789 2 0.6439 1 0.5292 1 0.4344 3 
24 0.4695 3 0.4849 3 0.4974 2 0.4483 1 
Alternative Ranking Using extended topsis 
mC] mC2 mC3 mC4 
aD 
en 
+ & 
© as nat 
—-ugsu e EA Soak 
Saco DN E $908 rats 
ko -3 i - = a aggeg 
Sgcs = Isgo 
T(1) T(2) T(3) T(4) 


Figure 2. Alternative Ranking using Extended TOPSIS Method 


Table 9. Final Scores of Alternatives 


TH TO TO TH 
cP 0.4661 0.4418 0.4806 0.4027 
c, 0.5117 0.5334 0.4967 0.4456 
Cc, 0.4789 0.6439 0.5292 0.4344 
c, (9 0.4695 0.4849 0.4974 0.4483 


5. DISCUSSION 


The problem of the selection of the head chef at a famous restaurant has 
been solved by applying the extended TOPSIS method based on the 
suggested distance and similarities measure. Based on the proposed method, 
four decision-makers were selected to select the best head chef to take 
charge of the kitchen. Each decision-maker assigned the rating values that 
corresponded to each attribute and alternative, as shown in tables 1-4. Using 
CIF-PIA and CIF-NIA, the separation measure between the alternatives was 


A. 
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computed, as shown in tables 5-6. Then, the last relative closeness 
coefficient was calculated, as shown in Table 7. Each value of the relative 
closeness coefficient was arranged in descending order and all alternatives 
were ranked corresponding to each decision-maker, as shown in Table 8 and 
Figure 2. According to T® decision-maker, Zə remains the PIA for the 
selection of head chef. According to T2? decision-maker, 33 is the PIA for 
the selection of head chef. According to TO? decision-maker, %3 is the PIA 
for the selection of head chef and according to T decision-maker, 34 is 
the PIA for the selection of head chef. However, if only one head chef needs 
to be selected, then the highest ranking PIA should be arranged in 
descending order and the best head chef should be chosen to take charge of 
the kitchen. 


5.1.Conclusion 


In this project, the issue of MCGDM under the CIFS-set environment 
has been discussed. An adaptation of the TOPSIS method has been 
illustrated to show the effectiveness of the proposed operators. According 
to the findings, these decision-making techniques can represent uncertainty 
more effectively than the current approaches and provide us with a 
comprehensive understanding of real-life scenarios. The problem of the 
selection of head chef for the newly opened restaurant has been dealt with 
in this research. The proposed approach of TOPSIS yields the best head 
chef. The findings of this research can be applied in the future to hypersoft 
set, interval-valued soft set, bi-polar soft set, and other ambiguous 
environments. 
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ABSTRACT 


Mathematical modeling is a vast field that has interdisciplinary implications 
for research. These models help to investigate the basic dynamics and 
quantitative behavior of infectious diseases that affect human beings, such 
as COVID-19, hepatitis B virus (HBV), and human immunodeficiency 
virus (HIV). The current study investigates the spread of HBV by using the 
basic virus model. In order to determine the stability of disease-free and 
endemic equilibria, the basic reproduction number is determined. The 
convergence and divergence of disease-free and endemic equilibria are 
demonstrated by using standard finite difference (SFD) and non-standard 
finite difference (NSFD) schemes. Arguably, SFD schemes, namely Euler 
and Runge-Kutta order four (RK-4) schemes, converge for lower step sizes, 
while the NSFD scheme converges for all step sizes. The latter is a strong, 
efficient, and reliable method that shows a clear picture of the continuous 
model. All the results are validated using numerical simulations in order to 
better comprehend the dynamics of the disease. The theoretical and 
numerical findings in this work can be applied as a useful tool for tracking 
the prevalence of HBV infectious disease. 


Keywords: convergence, divergence, HBV model, local and global 
stability, numerical schemes 


1. INTRODUCTION 


HBV is a devastating liver infection and an epidemic that poses a 
considerable health hazard, specifically in developing nations [1, 2]. HBV 
damages the liver and causes numerous shocking disorders. Almost 360 
million people worldwide are affected by this infectious disease and among 
them are more than 150 million HBV transporters (who transmit the disease 
to healthy people) [3, 4]. Also, more than 150 million of them have chronic 
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hepatitis C virus (HCV) infection and another 150—200 million are HBV 
chronic liver infection carriers [5]. HBV is transmitted by exposed sex, 
blood transfusion, and spreads into the recently born infant during 
pregnancy from the infected mother. 


Mathematical models [6-19] are very useful for analyzing the dynamics 
of real-world problems. Researchers utilize ordinary differential equations 
[6-11] as well as partial differential equations [12-15] to investigate 
different characteristics of physical objects. Fractional order [16-19] 
mathematical models are also among the most investigated topics in the 
world in this regard. Mathematical modeling also helps in the analysis of 
experimental results, as well as understanding the fundamental biological 
mechanisms of disease transmission. Many mathematical models have been 
suggested to improve the understanding and knowledge of HBV disease. 
Zeuzem et al. [6] and Nowak et al. [7] recommended simple virus infection 
models in which clear immune responses are not required. These models 
are mostly applied in the analysis of virus infection dynamics. In [8], the 
author proposed and examined the stochastic SACR model for the 
transmission of HBV disease. In [9, 10], the authors focused on population 
characteristics, infectious diseases, infection characteristics, and related 
social factors of HBV. Aja et al. [11] constructed a model based on 
statistical data which is not only useful for forecasting disease transmission 
but also helps to identify the key factors that affect the disease spread. 


In early 1980s, mathematical modelling was used to investigate the 
dynamics of HBV and the efficiency of disease control. To date, several 
systematic methods have been applied for the prevention of HBV 
transmission [20-22]. The most attractive and effective way of reducing the 
prevalence of HBV in newborns is vaccination. In case of vaccination, the 
emphasis shifts towards the use of a mathematical model to forecast the 
long-term consequences of immunization on hepatitis B regulation. Many 
scientists have analyzed the transmission mechanism of HBV and evaluated 
the effect of various vaccine policies [23, 24]. Recently, in [25], the author 
proposed and analyzed a non-linear basic adaptation system for HBV 
spread. The stability of disease-free and endemic equilibria was discussed 
for the continuous model. The author used a mathematical framework for 
the continuous model to study the stability of equilibria which strictly 
depend on the basic reproduction number. 
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The aim of the current study is to prevent the spread of HBV and explore 
its threats to the health of general people. Multiple numerical schemes, such 
as Euler, RK-4, and the NSFD scheme have been deployed to assess the 
various properties of the model and demonstrate their ecological feasibility 
and accuracy. The findings demonstrate that the aforementioned schemes 
give a thorough description of the continuous model. 


The structure of the current paper is organized as follows. The HBV 
disease model and corresponding parameters are described in Section 2. In 
Section 3, the expression for Rg, the equilibria of the system, and their 
stability are explored for the continuous model. SFD schemes, including 
Euler, RK-4, and the NSFD scheme are produced for the continuous model 
in Section 4. The fundamental dynamic properties of the continuous model 
containing the positivity of solutions cannot be preserved by Euler and RK- 
4 schemes for higher step sizes, leading to numerical solutions that are 
different from the solutions of the original system. However, the NSFD 
scheme gives positive solutions for all finite step sizes and basically 
overcomes the shortcomings of the SFD schemes. The analysis also 
demonstrates that the NSFD scheme is unconditionally convergent and 
provides better results in all aspects than Euler and RK-4 schemes. The 
theoretical findings are strengthened by numerical simulations at each stage. 
In Section 5, brief conclusions are offered to sum up the whole work. 


2. INTRODUCTION OF HEPATITIS B DISEASE MODEL 
Figure 1 depicts the non-linear basic HBV epidemic model [25]. 


P Ysw AF 


os pI uv 


Figure 1. Illustration of Hepatitis B Disease Model 


acp oS — ySV, 

dl 

a ySV — gl, (2.1) 
dV 

p = Al — pV 
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2.1. Parameters 

The parameters used in system (2.1) are described as follows: 
shows the number of uninfected cells in the population 
shows the number of infected cells in the population 
shows the number of free virus particles in the population 
shows the constant rate that explains the process reliability 
shows the constant rate 
shows the rate at which the virus leaves infected cells 
shows the die rate of uninfected cells 


shows how fast infected cells are dying 


* S&S 9 2» 0 cs5 wc à 


shows a free virus 
The parameters p, y, c, Q, A, u are all positive constants. 


3. STABILITY ANALYSIS OF CONTINUOUS MODEL (1) 
EQUILIBRIA 


The following subsection gives the equilibria of model (1). 
3.1. Equilibria of Model 


If all other classes are equal to zero aside from the susceptible class, the 
above described system (2.1) has the disease-free equilibrium (DFE) 
point Ey = (5, 0,0). If the system (2.1) is solved simultaneously, then the 
endemic equilibrium (EE) point occurs which is E* = (S, I, V*) = 
(GH, RD Ro- 

yA’ ya ? y 
the stability of the equilibria of model (1). 


3.2. Basic Reproductive Number (Rg) 


). In the following, Ro is provided which is crucial for 


The most essential threshold quantity for any infectious disease is Ro. It 
helps to assess whether a contagious disease will spread among the people 
or not. Rọ is a significant biological quantity that may be determined using 
the modified technique known as the “next generation matrix” [26]. The 
transmission matrix F (x) and translation matrix V (x) for system (2.1) can 
be described as follows: 
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FG) d maruse L , 2 " 


From the above, we can easily get 
ER 0 
ra c|anav = [*, | 
0 0 rs 
As Rọ = p(FV~?), so using the above information, Ry becomes 


— Yp^ 
Lupo” 
3.3. Stability of Equilibria 


0 


Theorem 3.1 [25]: The DFE point of system (2.1) is locally asymptotically 
stable (LAS) whenever Ry < 1, as displayed in Figure 2(a). 


Theorem 3.2 [25]: The EE point of system (1) is LAS whenever Rọ > 1, 
as displayed in Figure 2(b). 


(a) 


150 


100+ 


Population size 


50r 


0 20 40 60 80 100 0 100 200 300 400 500 
Time Time 


Figure 2. Numerical simulations for HBV model (1) with (a)y = 
.085, (b) y = .85. The remaining values of the parameters are offered as 
u = .196,0 = .0667,A = .001, 9 = .1, p = 10. 


According to the numerical simulations depicted in Figure 2(a), the DFE 
point of the continuous system (2.1) is stable if Rọ < 1. On the other hand, 
as illustrated in Figure 2(b), the EE point exists and becomes stable if Ro 
7]. It is evident from Figure 2 and the preceding discussion that the EE 
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point emerges and remains stable when the DFE point becomes unstable, 
that is, when Rọ passes the critical value Rọ = 1, the system exhibits 
transcritical bifurcation and stability exchange from one equilibrium point 
to another. It shows that Rp=1 is the saddle point for system (2.1). More 
details can be found in [25]. 


4. NUMERICAL ANALYSIS OF SYSTEM (1) 


In the current section, different discrete schemes for the continuous 
model (1) are constructed. The main aim is to better understand the 
dynamics of HBV disease transmission. Three schemes, namely Euler, RK- 
4, and NSFD schemes [27, 28] are constructed. It is demonstrated that the 
NSFD scheme is not only independent of step size but also more preferable 
as compared to the other defined schemes in every aspect. 


4.1. Euler Scheme 


Subsequent steps can be employed to construct the Euler scheme for 
system (2.1). 
spt. = SP 
h 
SP+1 — SP = h[p — oS? — ySPV?] 
SPtl = SP + h[ p — oS? — ySPV?]. 


= p — 0S? — ySPy? 


Similarly 
[p+1 — [v 
= yS?V? = gl? 
h 

JP*1 = [P + h[ySPV? — gI?], 
and 

+1_ 
a= AP — uy? , 


VP+ = y? + h[AI? — pV? J. 
Therefore, the Euler scheme becomes 
SP+ = SP + h| p — oS? — yS?V?] 
JP*1 = [P + h[ySPV? — gl] (4.1) 
VP+ = VP + h[AI? — uV?]. 
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Once the step size is small, then the numerical solutions offered by the 
Euler scheme are positive, as illustrated in Figure 3(a-c). According to 
Figure 3(d), stability 1s disrupted as the step size increases. Therefore, the 
stability and positivity of Euler scheme does not maintain for all finite step 
sizes. 


Population size 


Population size 


0 20 40 60 80 0 20 40 60 80 
Time Time 


Figure 3. Graphs generated by Euler scheme (4.1) with (a) h=0.001, (b) 
h=0.01, (c) h=0.1 and (d) h=3. The remaining values of the parameters are 
offered as y = 0.085, u = 0.196, ø = 0.0667,A = 0.001,09 = 0.1, p = 
10. 

4.2. Fourth-order Runge-Kutta Scheme (RK-4) 


To develop the RK-4 scheme, it is assumed that S = p,,] = q4 and V = 
rı. Then, the RK-4 scheme is described as follows: 
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Stage 1 

p, = hlp — (d+ yV?)S?], 
qı = h[yS"V? — gl], 

n = h[ AI? = nV? |. 


Stage 2 
sap (rev Ge) 
vip ere) -o( 2) 


n = nfa(?+ 2) -u (vo 2] 


Stage 4 
p, = h[p — (c + Y(V? 3)? + p3))], 
q4 = hly (V?  r3)(S? + p3) — eU + q3)], 
r, = h[AQP + q3) — u (V? + r3)]. 
Finally, we get 
1 
Set = SP + z [Pa t 2p; + 2p3 + pal, 
[P+ = - [q + 20; + 293 + q4], (4.2) 


ppt = In + 2rj + 273 + r,]. 
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(a) 


Population size 


Population size 


[2] N — Oo — Nw 
T T T T 


0 20 40. 60 $0 100 0 5 10 15 20 25 30 
Time Time 


Figure 4. Graphs generated by RK-4 scheme (4.2) with (a) h=0.01, (b) 
h=0.1, (c) h=1 and (d) h=3. The remaining values of the parameters are 
offered as y = 0.085, u = .196,0 = .0667, ^ = .001, 9 = 0.1, p = 10. 


Figure 4(a-d) provides a graphic representation of RK-4 scheme (4.2). 
The graphs produced by the RK-4 scheme are positive and convergent for 
small step sizes, as seen in Figure. 4(a-c). Figure 4(d) illustrates that the 
stability of model (1) is destroyed as the step size increases, that is, the RK- 
4 scheme diverges for large values of h. 


4.3. Non-standard Finite Difference (NSFD) Scheme 


Mickens [29] propagated the NSFD scheme which approximates the 
solutions of ordinary and partial differential equations. According to Shokri 
et al. [30], there are two factors that affect the study of the NSFD scheme. 
The first is how to discretize the derivative and the second constitutes the 
best way to approximate the nonlinear terms. The forward finite difference 
approximation is one of the general methods of first order derivative 


i mE >.. dy. ] : 
discretization. The first order derivative m is written in standard form as 


yGcrh)-y(x) 


: , where h stands for the step size. According to Mickens, this can 
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y(x+h)-y(x) 

9 (n) 
function. There are various formulas for (h) that can be viewed in [30]. 
In the current paper y(h) = h. The following steps can be taken to construct 
the NSFD scheme for system (2.1). 


be expressed as , where (h) is referred to as the denominator 


SP*1 — gp 
a =p- osPti — ySer ye, 
pet. = I? 
LE = ySPvyP = gree, 
p+1_yp 
——— = APP — yy? 
The above equations can be rewritten as 
1+ho+hyvP ? 

41 _ IP+hysPvP 
[P+ = EEUU (4.3) 
yeti = VP +hAIP 

~ jq4hp ^ 


4.3.1. Investigating the Stability of Equilibria for the NSFD Scheme. 
In order to discuss the local stability, suppose 


F = SPH! = ee 
1+ ha + heV?’ 
G = Ptt = ree (4.4) 
1-hqo 
— yp — VP +hAIP 
H=V E 


Theorem 4.1: The DFE point Ep for the discrete NSFD scheme (4.3) is 
LAS if Ro < 1. 


Proof: Consider the Jacobean matrix 


OF OF OF 
as OI OV 
0G OG OG 
_|ac ae ae) 4. 
J as OI OV i 
OH OH OH 
as OI OV 
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From (4.4), the following derivatives can be easily obtained 


OF _ 1 OF. OF . (hptS)hy 0G __ yw 0G | 1 0G . 

as  (i-ho-hyV)OI . '0V  [(i-ho-hyV)]2'8S —(1-hQ)O0I  (1-hg)8V - 
hyS OH 0H _ hA 8H 1 

(-hoy)8S ^8I  1+hu’ ôV  1+hp’ 


By replacing all the above derivatives in (4.5), we obtain 
1 


—— 0 0 
1-ho-chyV 
1 
J= 0 0 
1+hp 
hA 1 
0 
1thu 1+hu 


0 0 
1tho 
1 
J= 0 1+hy 0 
0 h^ 1 


1-hu 1+hp 


To verify the eigenvalues, we consider 


that 1s, 
ey | 0 0 
1itho 
1 
uy ^ " (4.6) 
0 hA 1 u 


1+hu 1-hu 


System (4.6) has the following characteristic equation 


(= a A) (x + z + =u) Pak) aim) = 0. 
(4.7) 


Using one of the eigenvalues from (4.7), A, = vus < 1 is obtained. The 
other two eigenvalues can be obtained from 
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M+ (s SR mu) A+ (1— Ro) Ferme 2 


According to Routh Hurwitz criterion [31, 32], the roots of the above 
equation are less than one if Rọ < 1. As a result, the DFE point for the 
NSFD scheme is LAS. 


Theorem 4.2: The EE point E* for the discrete NSFD scheme (4.3) is LAS, 
if Rọ > 1. 
Proof: In a manner analogous to Theorem 4.1, the Jacobian matrix is 
obtained as 

OF OF OF 


ðs ðL ðv 
aG aG ðG 
J=] or. wl (4.8) 
ðH OH OH 
ðs ðL ðv 


By putting all the derivatives from Theorem (4.1), the Jacobian matrix 
(4.8) becomes 


1 (hp+S)hy 
(1+ho+hyV) [((14+ho+hyv)]? 
_ yhv 1 hys 
L= (1+hg) (1+hg) (1+hg) Ve) 
hA 1 
0 
1-hu 1thu 
Furthermore, by replacing EE point E* in matrix (4.9), we obtain 
1 (hp+S*)hy 
(14+ho+hyV*) [(1+ho+hyv*)]? 
ee yhv* 1 hes* 
JE) = (1+hg) (1+hg) (1+hg) 
hA 1 
0 
1t+thu 1t+hu 
To check the eigenvalues, we take 
CE") — M] = 0, 
that is, 
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1 (hp+S*)hy 


(ho*hyV*) |. A 0 [ho RhyV*)]2 
yhV* 1 u hys* _ 
(1+hg) (1+hg) (1+hg) =0. (4.10) 
0 hA 1 
1+hu 1-hu 


From (4.10), A, = 


can be obtained from 
2_ 1 1 ) E ( (hp+S*)hy hA yhv* ) 
À "rr i 1+hu ( (1+ho+hyV*)? ^ 1thy ^ (1+hg) T 


( ictal Ya )) (Ry — 1) = 0. (4.11) 


1thu — (1+h@) 


1 . . " 
thot) 8 obtained. The other two eigenvalues 


The Routh Hurwitz criterion [31, 32] demonstrates that the roots of 
equation (4.11) are smaller than one if Rọ > 1. This demonstrates that E* 
for the discrete NSFD scheme (4.3) is LAS whenever Rọ > 1. 


(a) 


Population size 


Population size 


400 600 800 1000 05 2000 4000 6000 8000 10000 
Time Time 
Figure 5. Graphs generated by the NSFD scheme (4.3) with (a) h=0.1, (b) 
h=1, (c) h=100 and (d) h=500. The remaining values of the parameters are 
offered as y = .085, u = 0.196,0 = .0667, A = .001, o = 0.1 p=10. 
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Figure 5(a-d) depicts that the NSFD scheme (4.3) remains stable for 
each step size. The SFD scheme, however, shows convergence for lower 
step sizes. This demonstrates that for each step size the NSFD scheme is 
positive forever and unconditionally convergent. 


S. CONCLUSION 


The basic disease model for HBV has been proposed and examined in 
the current paper. For the proposed epidemic system, the equilibria and 
threshold quantity Rọ have been determined. In order to verify the 
convergence and divergence of the DFE and EE points, various numerical 
approaches have been employed. It has been demonstrated that RK-4 and 
Euler schemes depend on step sizes and are conditionally convergent which 
makes them ineffective. On the other hand, the NSFD approach developed 
for the current model is not step size dependent. This shows that the NSFD 
scheme is incredibly effective in comparison to the notable Euler and RK- 
4 numerical schemes. This scheme is a simple technique that demonstrates 
how discrete and continuous models behave appropriately and provide 
mathematically accurate results. By employing it, the spread of HBV 
epidemic diseases can be monitored effectively. The findings presented in 
this paper are beneficial to humanity in the field of medicine and can be 
applied also as a useful tool to predict the appearance of HBV epidemic 
disease. Each section includes numerical simulations to validate the 
theoretical findings. 
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Image Registration using the Rigid Group 

Muhammad Yousuf Tufail*, and Saima Gul 

NED University of Engineering and Technology, Karachi, Pakistan 
Abstract 


Image registration is the process of approximate matching of the source 
image to the target so that they resemble each other. In this study, two- 
dimensional image registration is presented using the rigid group. This group 
is a finite dimensional group (four-dimensional in this case) under 
composition. The dimensions of the rigid group are scaling, rotation, and 
translations along the axes. In this paper, an algorithm for the construction 
of rigid transformation is presented using the discretized objective function. 
This objective function is based on SSD (sum of the squares of the distances 
between the pixels intensities) and calculates the discrepancy between the 
images. The coarse search and the gradient descent approaches have been 
used for the optimization. The proposed algorithm is implemented on 
variety of images. The numerical examples illustrate the ability of the 
proposed algorithm. 


Keywords: algorithms, coarse search, image registration, optimization, 
rigid group. 


1. INTRODUCTION 


Image registration is the process of the approximate matching of the source 
image to the target so that their appearances resemble each other. The 
pioneer work of D’ Arcy Thompson remains the original inspiration in the 
field of image registration. In his famous book On Growth and form |1, 2], 
which was first published in 1917 and then re-published in 1942, he 
presented the idea of image transformation between species. Some of the 
images and their relative transformations that he referred in his book 
motivated the current authors to further explore finite dimensional group, 
which is the scope of this paper. Although there have been many 
developments in image registration since 1942, still the most remarkable 
and influential approach is to find a diffeomorphism between the two 
images using a method known as "Large Deformation Diffeomorphic 
Metric Mapping (LDDMM)" [3-7]. 
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Lisa Gottesfeld Brown published the first comprehensive review of image 
registration in 1992 [8]. Brown highlighted the various requirements of 
image registration namely (1) a feature space (contains information about the 
images), (ii) a search space (consists of a set of transformations), (iii) a search 
strategy (selection of transformations to get an optimal solution), and (iv) a 
similarity metric (to measure the discrepancy between the images). Brown 
suggested only three lists of transformations (all were finite) including (a) 
affine transformation (shearing, aspect ratio, and rigid), (b) perspective 
transformation (a transformation that maps 3D to 2D), and (c) polynomial 
transformation (suitable when minor information regarding camera geometry 
is available [2, 9]). The basic steps in image registration still remain the same. 
However, the choice of transformations has evolved drastically. The full set 
of diffeomorphisms, that is, non-rigid registration using infinite dimensional 
group is an integral part of image registration and widely applied in medical 
imaging [6—7, 10]. 


Several methods for image deformation have been developed and applied 
to the distinct fields of mathematics and computer science since 1992. These 
methods include artificial intelligence[11], image encryption [12], graph 
theory, [13] and many more [14—26]. The adoption of diffeomorphic image 
registration has caused a revolutionary and significant change in the field. 
Diffeomorphisms refer to those smooth functions whose inverse are also 
smooth [2, 9, 27]. 


In this paper, image registration is presented using the group of rigid 
diffeomorphisms. Two-dimensional grey scale images [2, 9] have been used to 
conduct the experiments. Image registration is subject to the optimal solution 
of the current study's objective function which is defined in Eq. (2.2). This 
equation demonstrates how close the images are. This objective function is 
based on SSD (sum of the squares of the distances between the pixels 
intensities). For more details, see equation (2.9) in Section 2.2 in [2]). 


2. RIGID GROUP AND IMAGE REGISTRATION 


The rigid group of diffeomorphisms is one of the most common group 
which is extensively used in image registration. It has been studied 
extensively over the past several decades, both for medical and non-medical 
applications [28-33]. This group consists of uniform scaling, rotations, and 
arbitrary translations along the axes [8, 19, 34—36]. Let us first present each 
of these separately. 
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2.1. Rotation 
Suppose an arbitrary point e in R?. Let ( € [0,27[ be the angle of 
rotation, then the rotation transformation w, to generate a new point (5) 
in IR? is defined as 
XN. XX. (cos € 2: 15 
" =a G) 7 s ( cos¢é b) 
2.2. Scaling 


Suppose y € R*,then the scaling transformation y; over È) is 


represented below as 


(P) -r() 
2.3. Translation along the Axes 


Suppose tr, tr, € R are the translations along the x-axis and y-axis 


x) is defined below 


respectively. The translation transformation p3 over S 


as 
X Xv (Xt tr 
ih TUM E) E (y + E 
The diffeomorphims group generated by the composition of scaling, 


rotation, and translations along the axes is termed as the rigid group. See 
Proposition 2.1 for more details. The composition of individual 


transformations yields (a) — uw (65) where Y = W3 ° Yz ° Wy 
Car(s OQ (E) ee anos en. 


The formula to represent the inverse rigid diffeomorphism is defined as 
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(y) = = (Sire rook (GE) = bale € [0,27[,0 # y € Rt 


H BOE) reete ren 


Proposition 2.1. The existence of the composition of rigid diffeomorphisms 
in the domain Q c R? implies the following fact: 


(1° Wr): (16) = pı: (p21), VU, V € rig(Q, R?). 

Here, '.' denotes the action of Y over the image I, o indicates the 
composition, see Section 1.2.4 in [2], and the set of rigid diffeomorphisms 
is represented with rig (Q, R?). 

Proof: Suppose 4 and w are rigid diffeomorphisms in Q such that their 
composition 4 o wp is also a diffeomorphism in Q, then 


(Wyo We) T(x) = A ov) Qo) 
= (z ° pD) 
= I(yz* (x) ° Yr (x) 
=] oz*(x) Y *(x) 
= (19v; Qoi Q9 
= (V; 160) opr) 
= yy, ° (V; - I(x)) 
Definition 2.1. Suppose J, andl, are two greyscale images. The 


discrepancy between the images is computed with the following matching 
functional or objective function 


MW) = f @ ey *)(x) - I(x) }’dxdy, x = (X,Y) ENC R^. (2.1) 


where y represents an element from the rigid group (See Section 2). Then, 
the process of finding an optimal w that minimizes Eq. (2.1) is termed as a 
rigid registration. It is clear that Eq. (2.1) is the continuous version of the 
objective function. As far as numerical computation is concerned, the 
continuous form is discretized. Suppose we define our discrete domain ID. 
This domain serves as the discrete domain for all the images presented in 
this paper. Suppose M is a positive integer and D = (([0, M — a]/(M — 
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a) — 0.5) x ([0, M — a]/(M — a) — 0.5); a = 1} isa uniform square grid 
that contains M? grid points, a simple example by setting M — 100 is shown 
in Figure 1. Consider M — 100 for all the numerical examples in this paper. 
The discrete form of the objective function (defined in Eq. (2.1)) is 


M M 
MG) 2 V Y (Gedy) - 0) - (2.2) 


i=1 j=1 
Here, yi; indicates the value of the image at position (i, j) in the image 
matrix. 


Uniform square grid 


0.5 


0.4 F 


0.3 F 


0.2 F 


0.2 0.4 0.6 


-0.6 -0.4 -0.2 
Figure 1. A Total of 10,000 Points for a Square Grid are Given. The range 
for Data Points is {[-B.n] X[-B.n];B =n = 3l with 100 Grid Points per 
Side, i.e., M — 100. 

Moreover, the values of /; are defined over the uniform grid points, 1.e., 
at yi; and not at vy (y j). Furthermore, it must be noticed that the domain 
of Ji » y! is not Q, instead it is y(Q) at the continuous level. It is highly 


likely (for some points) that J; » y ! move outside Q in the registration. 
Thus, Q is a restricted subset of R°. As far as computation is concerned, the 
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following values for Ji ° y (yjj) are used. 


1. If yw (y; j) EQ, bilinear interpolation (Section 2.1.2 in [2]) is used by 
using the values of /; corresponding to four neighboring edges of a square 
that contains y~! (y; j)- 


2. Ify (y; j) € Q,, a constant background value (0 for a black background 
in Example 4.1 and 1 for a white background in Examples 4.2 and 
4.3) is set. 


The central goal for the registration is to minimize M(v) for v €G. In this 
paper, two methods are used for optimization namely (1) the method of 
coarse search (that surveys the whole discretized domain for all the 
parameters, four in our case), and (ii) the trust- region-reflective method. 


The latter method is used when researchers encounter the non-linear least 
squares problem and it yields a local optimal solution. The trust-region- 
reflective algorithms have been used through MatLab, see example 4.3 . 


Eq. (2.2) denotes the general form of the least-squares optimisation 
(based on SSD) function for any choice of planar transformation group (1.e., 
G is a Lie group that acts on the plane [37], G is the rigid group for this 
study). This is a numerical optimization problem. Optimization is a huge 
field with many known algorithms [9, 38—40] whose applicability depends 
on the nature of the objective function. The purpose here is not to survey 
this field but to demonstrate some simple applications. 


3. COARSE SEARCH 


In this method, firstly, the domain for the matching functional is defined. 
This matching functional is computed for several points y € G, see 
Algorithm 4.1 for further details. This algorithm contains four nested 
loops, that is, the parameters of the rigid group. The values of individual 
parameters are set as Mz = 21, M, = 13, M, = 11,and M, = 11. The 
total search for optimal solution (minimization in this case) is the product 
of all four nested loops, that is, a total of 33,033 computations for the 
objective function. 


4. EXPERIMENTS 


In this section, the implementation of Algorithm 4.1 over synthetic data 
is presented. In synthetic data, the target is generated (by using known 
rigid transformation) from the source. Three examples are presented in which 
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the first example involves a pair of smooth images, whereas the other two 
examples are related to the registration of non-smooth images. 


Example 4.1 (pair of smooth images). In this first warm-up example, 2D 
Gaussian blob (an example of smooth image) is considered for image 
registration. The Gaussian function is considered as the source of image 
registration. A rigid transformation, with the parameters values 6 — 
pi/2, y = 0.8, tr, = 0.2 and try = 0.4, are used to generate the target. 
Algorithm 4.1 yields the optimized set of parameters whose values are 
Copt = T/2, Yopt = 0.8, txopt 7 02 and ty; = 0.4. These optimised 
parameters are used to transform the source with a rigid transformation. The 
transformed source is then subtracted from the target. By doing this, 
attention is drawn toward the fact that how good the images are aligned 
(transformed source and the target). Midgrey pixels indicate a perfect 
match,! whereas the evidence of mismatch appears as black (pixel value is 
= 0 ) or white (pixel value is ~ 1 ). Figure 2 displays the results of 
registration. An excellent registration is obtained (as expected). The optimal 
solution is highlighted with a red circle in the bottom part of Figure 2. 


o os ii 1 1.5 2 2.5 3 3.5 


Figure.2. Image Registration is Given for Example 4.1. Top: Left: Source 


! Perfect match here and elsewhere means good visual registration, i.e., there is no 
mismatch (apparently) between the transformed source and the target; and evidence of 
midgrey pixels screen when the transformed source is subtracted from its corresponding 


target. 
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(the Gaussian exp (—5x? — 3y?)), Target, Deformed Source, and 
Difference Image (Gh oy 1—1L-41)/ 2), Uniform Mid Grey Confirms 
that Perfect Registration between the Images is Obtained; Right: Rigid 
Diffeomorphism y~* Applied to the Target Grid with the Comparison of 
Contour Plots of I, o ij ^! (the Transformed Source ) and I, (the Target). 
The contour plots Have Been Made by Plotting Green Contours on Top of 
Red Contours. No Red Line Can Be Seen Due to Perfect Registration. 
Bottom: The Registration Errors Are Given. The Optimal Solution Is 
Identified as A Red Circle. 


Example 4.2 (pair of non-smooth images). Image registration for non- 
smooth images is always challenging [9]. Consider the pair of non-smooth 
images. Alike the previous example, again the synthetic data is considered, 
that is, 6 = pi/2,y = 1.1, tr, = —0.2, and try = 0.1 are set as the initial 
parameters for the known rigid transformation. The implementation of 
Algorithm 4.1 produces optimized parameters that are Coy = = Yopt = 


1.1,t = —0.2, and t — 0.1. After obtaining these optimized 


Xopt Yopt 


parameters, the source is transformed by following the same steps as in the 
Example 4.1. The results of image registration are given in Figure 3. Alike 
the previous example, a perfect registration is obtained. 


source 


ie 
RERE 


HH 


: | ! 
Hu 


i" 


t 
f 
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7 


Figure 3. Implementation of Algorithm 4.1 Over A Pair of Nonsmooth 
Images Used for Example 4.2. 
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Example 4.3 (pair of non-smooth images). Now consider a slightly harder 
example in which the target is generated from the source with the help of 
rigid transformation. The values ¢ = 1.6, y = 1.25, tr x = —0.27 and tr, = 


0.13 are set as parameter values for rigid transformation. Resultantly, 
Algorithm 4.1 produced Çopt = 1/2, Yopt = 1-2, txon = —0.3 and t 


opt Yopt 
0.1 as optimized parameters. These parameter values are noticeably 
different from those values of parameters that were used to construct the 
target. The reason for discrepancy (in values) is due to the fact that the set 
of allowable values for the parameters in coarse search do not contain the 
exact values of the parameters that were used to generate the target. Thus, 
Figure 4 displays adequate but not perfect registration. 
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Figure 4. Result of Image Registration Using Algorithm 4.1 for 
Example 4.3. 
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MatLab based optimizer Isqnonlin is used for the current non-linear 
problem. This optimizer is used for the (typical) least squares functions for 
optimization [2, 41]. 


Isqnonlin requires an initial guess for optimization. For this purpose, the 
above optimized parameters that were obtained at the first step of 
optimization (the coarse search), that is, Çini = 1/2, Yini = 1.2, txipi = 
—0.3, and ty, , = 0.1 were chosen as an initial guess. The optimizer yields 
a global minimum and the corresponding values of the parameters are 
Gopt = 1-6, Yop = 1.25, ty, = —0.27, and ty ope = 0.13. Figure 5 confirms 


t t 


a perfect registration. 
Original Source Target ,  Tunstormed grid , Target arid 
6 
> 2h, - 0.2 yy 
o 
02 ji 0.2 
Transformed source Fixing the Source over Target 


Figure 5. Results of Image Registration for Example 4. 3 Using Lsqnonlin. 
Top Row: Pair of Relevant Images for the Registration. Second Row: The 
Transformed Source (1,°y ! ) And the Subtraction Between the 
Transformed Source and The Target. Bottom Row: The Transformed Grid 
and The Target Grid. These Identical Grids Indicate the Perfect Registration 
Between the Source and The Target. 


Algorithm 4.1. Minimization of Eq. (2.2) using Coarse Search 
Input : Source J, and target I; 

output: Warp y ! along with the transformed image I, » iy! 
for ¢ = 0: 7/10: 27 do 

for y = 0.1:0.1: 1.3 do 
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for tr, = —0.5: 0.1: 0.5 do 

for tr, == —0.5:0.1:0.5 do 

Computation of deformation of source using bi-linear interpolation, 
loy !(xj)vyi ES 

Compute d = ||I, » ij^ (xij) — B(x) Yyy ES 

Compute y^! = i^! for the optimal solution of d, 

final deformation of source: I, o ^! 


5. Conclusion 


In this research, an algorithm for the construction of rigid 
diffeomorphisms has been presented. This algorithm gives an optimal 
solution for the objective function which calculates the discrepancy 
between the images (see Eq. (2.2)). It worked very well for the selected 
examples and found the optimal solution for perfect match. However, 
there is a deficiency in the algorithm that it gets stuck on the local 
minimum for some examples (not reported here). Due to this deficiency, 
the current authors were unable to find a perfect match (see Example 
4.3). Thus, other optimization techniques are also needed to overcome 
this deficiency. Gradient descent (there are few more, see [42, 43]) is 
one of the approaches and used in Example 4.3. 


The major benefit of the proposed algorithm is that (as it is based on 
coarse search) it surveys the whole domain and provides a good initial 
guess for more sophisticated optimization. Apparently, this method is 
very simple but it is always recommended to search for each parameter 
(over a relatively small number of values)in the whole domain. In light 
of modern research on optimization, it is worthy to implement coarse 
search because (1) it 1s worthless to compute the derivative at the first 
step (function might not be differentiable) and (ii) function might have 
several local minima scattered all over the domain. It becomes nearly 
impossible for the derivative-based optimization method (such as 
steepest descent) to find the global minimum. In this situation, the 
method of coarse search is the best option to get a better initial guess 
(possibly very near to the global minimum, see Example 4.3) in order 
to find the desired optimal solution for the perfect match. 
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ABSTRACT 


This paper is concerned with the numerical methods of susceptible exposed 
infectious recovered (SEIR) epidemic model of coronavirus disease 2019 
(COVID-19). The model is explicated numerically with three numerical 
schemes, forward Euler, Runge-Kutta of order 4 (RK-4), and the proposed 
non-standard finite difference (NSFD) technique, respectively. In the 
epidemic model of infectious diseases, positivity is the main property of a 
consistent framework, since the negative value of a subpopulation is 
useless. The NSFD technique ends up being a more important and trustable 
numerical system than forward Euler and RK-4 techniques since it 
preserves positivity, stability, and convergence. On the contrary, forward 
Euler and RK-4 schemes do not hold these characteristics for some choices 
of step sizes. Numerical simulations confirmed the findings. 


Keywords: convergence, NSFD method, SEIR model, stability 
1. INTRODUCTION 


In 1965, scientists identified the first human coronavirus case [1]. Later 
on, experts discovered a group of related human and animal illnesses caused 
by a family of viruses named coronaviruses after their crest-like appearance. 
People can contract seven types of coronavirus infections. COVID-19 
sickness is a serious illness brought on by the SARS-CoV-2 infection [2]. 
The SARS-causing strain first appeared in southern China in 2002 and 
quickly spread to 29 different nations. By July 2003, 774 individuals had 
died and more than 8,000 had been infected [3]. Fever, headaches, and 
breathing problems, such as coughing and shortness of breath, are all 
brought on by coronavirus infection [4]. 


Middle East Respiratory Syndrome (MERS) first appeared in Saudi 
Arabia. Many people who were exposed to it developed difficulty in 
breathing. It was difficult to treat the respiratory illnesses of such people 
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and they could not recover without needing special care. People who have 
hidden illnesses, such as diabetes, cardiovascular disease, or risky 
development, are more likely to transfer recognized illnesses. Being 
educated about the condition and how the disease spreads is the greatest 
strategy to prevent and tone down its transmission. 


On December 31, 2019, WHO received information regarding instances 
of pneumonia in Wuhan City, China with enigmatic explanations. On 
January 7, 2020, Chinese experts identified an original coronavirus as the 
cause and gave it the temporary designation ‘2019-nCoV’ [5]. According to 
[6], USA and Brazil are the two nations that are most affected by COVID- 
19. Dhandapani et. Al provided stiff, fuzzy IRD-14 day average 
transmission model of COVID-19 pandemic disease [7]. To examine the 
factors involved in the spread of the coronavirus pandemic in Pakistan and 
its potential controls, S. Ullah and M.A. Khan developed a numerical model 
[8]. Duccio and Fanelli focused on the ephemeral aspects of the coronavirus 
disease 2019 in China, Italy, and France for the time period January 22 - 
March 15, 2020 [9]. The focus of Xiao and Ruan [10] remained on a 
pandemic model without a monotonic occurrence rate. The prescribed 
model depicts the psychological impact of specific illnesses on the local 
population as the number of infectious agents increase. 


Since the spread of coronavirus has resulted in a pandemic, numerous 
mathematicians have conducted various analyses to create a model to 
predict its expansion [11—29]. 


The NSFD scheme was introduced by Mickens [23]. The construction 
of the NSFD mathematical model for a two-layered differential condition 
was the focus of Cresson and Pierret [24], who also investigated the model's 
various aspects, such as combination and solidness of the plan, among 
others. Selected mathematical models were established and tested using the 
RK method. Euler algorithm for requests 2 and 4 were completed. In order 
to account for a jungle fever model, Anguelov et al. [25] constructed the 
NSFD plot and investigated its dependability and constituent parts. To 
resolve the nonlinear Riccati differential condition, Riaz et al. [26] 
presented a very stable non-standard limited distinction (NSFD) strategy. 
By comparing the results to those of other mathematical operations, such as 
Euler and RK-4, the proposed model's accuracy and productivity were 
verified. A dynamic model for infection transmission components was 
proposed by Rafiq et al. [27]. 
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A mathematical model is developed in the current study by comparing 
the results with prominent mathematical schemes including the Runge- 
Kutta approach and the Euler scheme. In order to investigate how insect 
vectors affect the development of plant diseases, Rafiq et al. [28] proposed 
a one-of-a-kind, genuinely stable, non-standard finite differentiation 
(NSFD) numerical arrangement. Numerical assessments were supplied and 
the outcomes were distinguished. In order to plan the proposed model and 
determine the value of the limit boundary Rọ for the model, Rafiq et al. [29] 
developed a model consisting of nonlinear differential conditions. To 
analyze the susceptible exposed infected recovered (SEIR) dynamic model 
of COVID-19 (coronavirus) with reference to bimodal infection 
transmission in vulnerable populations, Ahmed et al. [30] proposed a 
construction protecting nonstandard limited distinction plan. The suggested 
mathematical framework produced workable arrangements for the 
perplexing bi-modular SITR nonlinear model, combined actually consistent 
states, and reflected dynamic consistency with the perpetual feel of the 
model. The model's analysis revealed that it maintained some stability at 
focused consistent states. For the mathematical setup of the SEIR response 
dispersion pandemic model, [31] provided two innovative NSFD schemes. 
Moreover, [32] proposed an SIR model in a fuzzy environment. Euler, 
Runge Kutta of order 4 (RK-4, and the NSFD methods respectively were 
developed with fuzzy extensions for the solution of the model. Dayan et al. 
presented rumor based fuzzy model and developed an NSFD scheme for its 
solution [33]. Some researchers studied fuzzy epidemic models using NSFD 
schemes [34, 35]. Ali et al. applied the NSFD scheme for the numerical 
solution of a cancer disease model [36]. Alsallami et al. used Euler, RK, and 
NSFD methods to solve a real-world problem [37]. 


The remainder of this article is structured as follows. Model formulation 
is presented in Section 2. Section 3 discusses the model's equilibrium states, 
stability, and threshold analysis. Section 4 presents numerical simulation 
and discussion of the results, while Section 5 displays the numerical results. 
The novelty of the current work lies in the fact that the proposed model has 
not been analyzed before using the NSFD method. 


2. SEIR MODEL FOR COVID-19 


Consider the model 
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dS BSI 

rx A—yS- WT (2.1) 
dE BSI 

ae a (ut+e)E (2.2) 
dl 

an — (y* ua) (2.3) 
dR 

PUE yl — uR (2.4) 
2.1. Theorem 


The system (1 — 4) has a DFE point (So, Eo Jo Ro,) = G ,0,0,0) 


BSI 
- ENS 2 
m '(ure) ’ 


and an endemic equilibrium point (S*, E*,1*, R*) = ( 


E I . 
E z) respectively. 


y+u+a u 
2.2. Basic Reproductive Number 


By incorporating the next generation matrix (NGM) approach, the value 
of Ro is determined. Let X = [S, E]* , then n = y(x) — z(x), where 


BSI 
v 
0 


and 


E (u + e)E 
20) = m +u +a) - E 


Y and Z are now the Jacobeans of the Y (x) and z(x) respectively at 
DFEP. They are listed below. 


B ^ 
C8) 
0 0 


and 


School of Science ie 
Volume 7 Issue 1, 2023 Fe 


Construction and Analysis of a Nonstandard... 


nee y+ i + x 


Inverse of matrix Z is 


: 0 
Z1 = (u F €) 
€ 1 


(ut+el(ytuta) (y*tpca) 


Computing the product Y Z7! we get, 


1 
B^ — 0 
eb [0 oP (u + €) 
0 0 


Qr ey rae) (Pare) 


Be B 
YZ = (eco ee) 


0 0 
Calculating the spectral radius of YZ~* we get, 
€ 
1 p 


"(yt otia) 
The spectral radius of YZ ! is equable to Rọ which is the maximal 
eigenvalue of YZ 1. 


Therefore, 

= Be 

(t 9otuta) 

2.3. Normalized Forward Sensitivity Index (NFSI) 


(2.5) 


Ro 


Chitins developed the idea of sensitivity indices by calculating the 
sensitivity indices of Rg [38]. A parameter's normalized forward sensitivity 


index (NFSI) is calculated as NFSI(¢) = RE (Ro). 


The NSFIs of f, €, B, y, and a are clalculated as follows: 


BORo _ B a(Be) MEDICI € E 
Ro0B —__BE___ (yu e) tria) € (u+ €)(ytu+@) ' 


(u +E€)(y+u+a) 
Vz 
92 — SIR 


NFSI(f) = 
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Similarly, 
NES) = — gp — süproHed = G+ 
(u + €)(y+p+a) 
NFSI(u) = — ENS ERECT, and 
NFSI(a) = eum 


These findings unambiguously demonstrate that f is the most sensitive 
parameter. 


2.4. Stability Analysis 
The Jacobean of the system (1-4) is 


d 0 2 0 
u-By Px 
E I S 
J=| By -&*9 p 0 
0 —€ —(y+u+a) 0 
0 0 Y =H 
u 0 —p 0 
I 
o= r cem wi : 
0 —€ —(y+u+a) 0 
0 0 Y =H 


If the absolute eigenvalues of J (o) are smaller than unity, that is, |A; | 
< 1,i= 1, 2,3. The numerical scheme converges in all cases. The eigenvalue 
from the aforementioned Jacobean matrix were extracted, where A, = A47 
—p, A, 2—(p + €), and 452 —(y + p + a). 


2.4.1. Stability at Endemic Equilibrium Point. Jacobean matrix at 
the EE point is 


0 zi 0 
u-Py dr 
jen=| gl ase = J 
By ope Pw , 
0 —€ —(y*tu-*a) 0 
0 0 Y -j 
where 
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S" = 5 and I’ E 
ut pu (y*u*a) 


Using the MATLAB database, the biggest eigenvalues of /(E*) were 
plotted. At the endemic equilibrium point, the Jacobean spectral radius has 
a maximum value that is smaller than unity, as is evident from the preceding 
figure, proving the intended assertion. 


3. NUMERICAL MODELING OF SEIR MODEL FOR COVID-19 
3.1. Forward Euler Method 


pI” 
Ce S” + h[A — uS” — S" =~] (3.1) 
S" 

Eu =E” +h É — (u+ €)E” (3.2) 
T^t12]" +h [EE* — (y + u + o)1" ] (3.3) 
R™+1 = R” +h [yI” — uR”] (3.4) 

3.2. Fourth Order Runge-Kutta (RK-4) Scheme 

3.2.1. Step 1 
pI” 
kı = S” + h[A — ps? — S™——] (3.5) 
pa 
m, = E" + hl —(ut+e)E”] (3.6) 
n, —I"-4h[cE" —-YYrpta@i"] (3.7) 
pı = R” +h [yI" — 4 R"] (3.8) 
3.2.2. Step 2 
k n 
kk OUS WS 
kg = h[A— gu (S" + -> — B ————————— (3.9) 
2 N 
(sn + ym e BY - 

mz = h[B——— p - # EE * 7] (3.10) 

m, T4 
n; = h [e(E” T3 — (y + p a)" uic! (3.11) 
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n D 
p; = h [y0 € 7) - u(8" 4 7)] 
3.2.3. Step 3 
"uu Fay (n + 
k= hasn er ——— 
k n 
(S7 4 79) +5) 


ms = h[B m 


m»; 
- Qt ec" +] 
m, Ny 
n3 —h [(E" + 7) — (y seit ace” *3)] 


E ng 2) n ,D2 
p= h [rne T) m 

3.2.4. Step 4 

S" + ky" + 
ka = h[A - a(S” + ks) -gg tX tr 
(S? + k3)(I" n3) 
N 

n4 =h [e(E" * m3) — (y + u + a)” +73)] 
p, = h [yU" t n3) — pCR" + p3)] 

3.2.5. Step 5 


m, = A[B — (u + €)(E” + m3)] 


1 

SH $c [ki 2k; 2ks + ke] 
1 

Eee = E” +z [my + 2m; + 2m; € m, | 
1 

pet IN +o [ny + 2n; + 2n; + n4] 


1 
RM = OR" += [pi + 2p2 + 2p3 + Pal 


3.3. Non-Standard Finite Difference (NSFD) Scheme 


(3.12) 


(3.13) 


(3.14) 
(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


To develop an explicit NSFD scheme, consider the above system as 
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dft) f(t + h) — f(t) 
dS — h 


S”, E". I”, and R "are approximations. Here, h is the time step size. 


+ 0(R) h Q0. 


gnti = gn BI" 
Se n+1 — gntin_ 3.25 
h uS N ( 3.25) 
os S” hA 
$a = a (3.26) 
lt hut h^- 
fe ps" 
Enti LO ME (3.27) 
14 h(u 4 €) 
pei I" +h eg? os 
~ 1+h(y+u+a) (pn) 
R” +hyI” 
R™*1 = ————— 3.29 
1+hu ( ) 


3.4. Stability Analysis of the NSFD Technique 


The NSFD scheme of the SEIR model is carried out in this section at 
DFEP (So, Eo Ig Ro) = G , 0, 0, 0). Consider the Eq. (3.25) to (3.28), 


m ntl — S+hA 

F=S Bu (3.30) 

- E+ nr 
G=E = — 3.31 
14 h(u 4 €) m 

I+heE 
H =p" = ME 3.32 
1+h(ytuta) oe 

R+hyl 
K = pt = m (3.33) 


Jacobean matrix of the Eq. (3.30) to (3.33) at the DFEP is as follows: 
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1 GERA) 
14h 7 [1+hp]? 
0 —— apa 0 
J(£o) = 1+h(u+e) UN 
he 1 
1+h(ytuta) = 1+h(yt+pt+a) 
0 0 hy 1 
1+hu 1+hp 


The eigenvalues of the aforementioned Jacobean matrix are 4, = ius < 


1 and Az = ez < 1. The following matrix is used to compute the remaining 


eigenvalues. 
1 hpA 
| 1+hu+e) uN 
JGo) = he 1 


1+hytuta) 1+hytuta) 


The other eigenvalues are quite complicated to be calculated algebraically. 
Therefore, they have been calculated numerically (shown in Figure 1). 


Disease Free Equilibrium 
aa -y -— -— T 


O.89 f T " i 


T E ] 


t 
So 


0.88 


0.87 


0.86 


0.85 f 


0.84 


Spectral Radius of Jacobean Matrix J(£ ) 


s " E 
o 10 20 30 40 50 60 7O 80 90 100 
Step Size 


Figure 1. Eigenvalues of J (čo) 


The spectral radius of the Jacobean matrix at the DFE point has a 
maximum value that is less than unity, as is evident from the following 
figure, proving the intended assertion. 
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4. NUMERICAL SIMULATION OF THE SEIR MODEL 


The findings are discussed in detail in this section. Additionally, Euler, RK- 
4, and NSFD methods in the model are compared in terms of their use. 
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Infected population 
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Figure 2. Infected Population using Forward Euler Scheme at h =0.0 
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Figure 3. Infected Population using Euler Scheme at h =0.1 


Scientific Inquiry and Review 


5 — SIR 


Volume 7 Issue 1, 2023 


Dayan and Iqbal 


we 410304 Disease Free Equilibrium 


10 


^ o % 


Infected population 


N 


0 5 10 15 20 25 30 35 40 45 50 
Time (days), h=1 
Figure 4. Infected Population using Euler Scheme at h =1 
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Figure 6. Infected Population using RK-4 Method at h =0.01 
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Figure 7. Infected Population using RK-4 Method at h =0.1 
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Figure 8. Infected Population using RK-4 Method at h =1 
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Figure 9. Infected Population using RK-4 Method at h 210 
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Figure 10. Infected Population using NSFD Scheme at h =0.01 
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Figure 11. Infected Population using NSFD Scheme at h =0.1 
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Figure 12. Infected Population using NSFD Scheme at h =1 
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Figure 13. Infected Population using NSFD Scheme at h =10 
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Figure 14. Comparison of Euler, RK-4, and NSFD Scheme at h=0.01 
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Figure 15. Comparison of Euler, RK-4, and NSFD Scheme at h=0.1 
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Figure 16. Comparison of Euler, RK-4, and NSFD Scheme at h=1 
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Figure 17. Comparison of Euler, RK-4, and NSFD Scheme at h=10 
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The behavior of the SEIR model is presented in the above graphs. Figure 
2 and Figure 3 show that the Euler technique behaves favorably and converges 
at small step sizes (h = 0.01 and h = 0.1). In Figure 4 and Figure 5, the 
Euler method shows divergence at slightly large step sizes (h = 1 and h = 
10). In Figure 6 and Figure 7, the RK-4 method shows convergence at small 
step sizes (h = 0.01 and h = 0.1). The method shows divergence as the 
step size increases, as shown in Figure 8 and Figure 9, respectively. In Figure 
10-13, the NSFD method converges to the same equilibrium points at 
different step sizes. In Figure 14-17, a comparison of Euler's method, RK- 
4 method, and NSFD method is depicted. It is clear from the above graphs 
that all models show similar behaviors at small step sizes and converge to 
the same equilibrium points. The Euler and RK-4 methods show divergence 
and negative behavior as step size increases, while the NSFD method 
converges to the same point. Graphical depiction of behaviors demonstrates 
that the Euler and RK-4 methods only provide convergence solutions for 
small step sizes and remain inconclusive for large step sizes. The NSFD 
scheme, on the other hand, exhibits good behavior and also provides a 
convergent solution for step sizes that are extremely large. 


5. CONCLUSION 


The mathematical study of the SEIR model for the spread of COVID- 
19 was carried out in this research. For this purpose, the reproduction 
number Ro, sensitivity of Ro, and equilibrium points of the model were 
determined. It was shown that both equilibrium points have similar stability 
properties. In order to solve the investigated model, Euler, RK-4, and NSFD 
algorithms were deployed. Different time step sizes were used for numerical 
experiments at DFE and EE locations. The collected results were examined 
and compared. From the said examination, the current study concludes that 
the suggested method yields findings that converge to true stable states for 
any time step size. However, the Euler and RK-4 methods do not hold for 
large time step sizes. Moreover, the NSFD method is bounded, dynamically 
consistent, and preserves the positivity of the solution which are important 
requirements when modeling a prevalent disease. 
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